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I. INTRODUCTION 

The studies of black-hole binaries (BHBs) that imme- 
diately followed the 2005 breakthroughs in numerical rel- 
ativity [IH3] soon revealed the importance of spin to the 
orbital dynamics [3]. One of the the most striking re- 
sult was the unexpectedly large recoil velocity imparted 
to the remnant due to an intense burst of gravitational 
radiation around merger [5J [5] . Recoil velocities as large 
as 4000 km s -1 were predicted for maximally spinning 
black holes [7] (in a configuration with both spins lying in 
the orbital plane, known as the superkick configuration). 
This prediction was based on model for the recoil veloci- 
ties that included the terms linear in the in the individual 
spins of the merging holes [5j |8] . These predictions trig- 
gered several astronomical searches for recoiling super- 
massive black holes as the byproduct of galaxy collisions, 
producing several dozen potential candidates [9HI4]. 

Accretion effects p~5j [16] would tend to align the spins 
of the BHs with the orbital angular momentum, sup- 
pressing the superkick and, apparently, the likelihood of 
observing large recoils. We recently found [TTJ HH] how- 
ever, that there are nonlinear spin couplings that lead to 
even larger recoil velocities when the spins are partially 
aligned with the orbital angular momentum. These so 
called hangup kick recoils can be as large as 5000 km s _1 
(see Fig.[l]). 

In this paper we continue our exploration of unex- 
pectedly large nonlinear contributions to the total re- 
coil [TTJ HH] . Here we concentrate on equal-mass black- 
hole binaries (BHBs) that precess. Our ultimate goal is 
to derive an empirical formula that takes into account 
all contributions to the recoil (at least to the level of a 
few percent accuracy). This would be a near hopeless 
task if we just start with a set of random configurations. 
Rather, we propose a program for developing sets of con- 
figurations with exact or approximate symmetries that 
allow us to model the recoil term by term. For example, 
in the hangup kick configurations |17[ 118) , the BHBs can 
be described by two parameters, the z component of the 




FIG. 1: The hangup kick configuration. Here Si z = S^z, while 
Six = —S2x and S\ y — —S2y The hangup kick configurations 
are preserved exactly by numerical evolutions. 



total spin S z , and the in-plane component of A (S2 — Si 
in the equal-mass case). If the generic recoil is also a 
function of A 2 , those terms would be suppressed in the 
hangup kick configurations. Here we continue the explo- 
ration by evolving configurations that activate different 
possible terms for the recoil. 

Not all nonlinear terms in the spins lead to large in- 
creases in the recoil [T^l [2D] . We therefore need to per- 
form many diverse simulations to try to elucidate which 
nonlinear terms contribute significantly and which can 
still be ignored. Here we explore the effects of preces- 
sion on recoils. We perform simulations of equal mass, 
precessing binary black holes, but also discuss the more 
general case of unequal masses. We also extend the phe- 
nomenological formulas for predicting recoils to include 
higher powers of the spins, explicitly including up to 
fourth order, making use of discrete symmetry proper- 
ties of the BHBs. 
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II. HIGHER ORDER RECOIL VELOCITY 
EXPANSIONS 

Our approach to the modeling of recoil velocities of 
merged black holes is based on the numerical evidence 
that the vast majority of the recoil is produced by the 
anisotropic emission of gravitational radiation at the very 
last stage of the merger, i.e., when a common horizon 
forms (see, e.g., Fig [3] below) . Around that time of the 
merger a burst of radiation of short duration occurs when 
spinning black holes merge [21] ■ In order to model the 
dependence of the recoil on the spins of the individual 
holes (and the BHB's mass ratio) we were guided by the 
leading-order post-Newtonian (PN) expressions for the 
instantaneous linear momentum radiated |22j (higher- 
order PN couplings can be found in [23]). Although, 
clearly the PN approximation is not valid at the moment 
of the merger (PN theory does not even account for hori- 
zons), our ansatz is that the parameter dependence in 
these PN expressions yields a useful starting point for 
constructing empirical formulas for the recoil. Thus a 
PN expression of the form F(f, P)n- S, where n is a unit 
vector in the orbital plane, becomes a fitting term on 
our formula of the form AS± cos(4> — fa), where A and 
fa are fitting constants. We then verify these assump- 
tions a posteriori by comparing our predictions to recoil 
results from other simulations. This is the basis of our 
phenomenological approach to modeling of recoil veloci- 
ties and can be summarized in an expression for the three 
components of the velocity linear in the individual spins 
of the holes [3 HUE] 

Recoil (q, a) = v m ex + v±_ (cos £ e\ + sin £ e 2 ) + v\\ h\\ , 
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and r\ = q/(l + q) 2 , with q = m\/m 2 the mass ra- 
tio of the smaller to larger mass hole, <3; = Si/m 2 , 
rrii is shorthand for mm the mass of BH i, the index 
_L and || refer to perpendicular and parallel to the or- 



bital angular momentum respectively, ei,<32 are orthog- 
onal unit vectors in the orbital plane, and £ measures 
the angle between the unequal mass and spin contribu- 
tion to the recoil velocity in the orbital plane. The an- 
gles fa\ and <ps are defined as the angle between the 
in-plane component A 1 - = M{S^ jm 2 — Sf~/mi) and 
S = + S 2 respectively and a fiducial direction at 
merger (see Ref. [25] for a description of the technique). 
Note that A = M(S 2 /m 2 — B\/m\) can be expressed as 
A = M 2 (c? 2 — qdx) I (l+q)- Phases fa and fa depend on 
the initial separation of the holes for quasicircular orbits 
(astrophysically realistic evolutions of comparable masses 
BHs lead to nearly zero eccentricity mergers). 

The current estimates for the above parameters are |251 
[26] : A m = 1.2 x 10 4 km s" 1 , B rn = -0.93, H = (6.9 ± 
0.5) x 10 3 km s" 1 , K = (5.9 ± 0.1) x 10 4 km s _1 , and 
£ - 145°, and K s = -4.254. Here we set B H and B K to 
zero, which is consistent with the error estimates in [24] 

Additional corrections from the hangup kick and effect, 
and a new effect, which we will dub the cross kick effect, 
are examined here. We note that these new effects were 
found using equal-mass BHBs, thus their dependence on 
mass ratio is still speculative. 



A. Symmetry considerations 

In our recent studies of BHB mergers we have found 
that nonlinear terms in the spin play an important role 
in modeling recoil velocities [17] . Here we will inves- 
tigate higher-order models for the recoil velocity based 
on the symmetry properties of its components [27] . In 
Boyle et al |27j a new method for developing empirical 
formulas for the remnant BH properties was proposed. 
This new method was based on a Taylor expansion us- 
ing symmetry properties to limit the total number of 
terms. Here, we will combine the two methods by us- 
ing PN-inspired variables for a Boyle et al type of ex- 
pansion. Fundamental to this construction is the be- 
havior of the BHB under discrete operations such as 
exchange (X) of the black holes labels (1 < — > 2) and 
parity (P) (x — > —x,y — > —y,z — > —z). Here we will 
model higher-order contributions to the recoil using the 
variables A = M(S 2 /m 2 - Si /mi), S = Si + S 2 , and 
i] = q/(l + q) 2 (where q = 7711/7712), for a closer connec- 
tion with the PN expansions [23] and for continuity with 
our previous parametrization as in Eq. ([2| above. 

Note that we could use a "simpler" set of variables, 
such as S±/m 2 = (Si±S2)/m 2 and 6m = {m\ — m 2 ) /m, 
with 777 = 777! + m 2 . We could then re-absorb the differ- 
ences in the definition of variables by allowing coefficients 
in front of each expansion term to depend on the mass 
ratio. For instance, in the V m case 



and because 



V m = a. 5m + b.Sm + 



1] = (1 - <5tt7 2 )/4, 



(3) 



(4) 
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TABLE I: Symmetry properties of key quantities 



Symmetry 




X 


o±/m =(bi + b2)±/m 






o\\/m = (bi + b2)\\/m 


+ 


+ 


A±/m 2 — (S2/m>2 — Si/mi)±/m 


+ 




Ay/m 2 = (S 2 /m 2 - Si/mi)\\/m 




+ 


n = fx — ?2 


+ 




5m — (mi — rri2)/m 


+ 




Vi_ 


+ 




n 




+ 



this is equivalent to the more usual (See Eq. ^ above) 
V m = if.5m.{A + B.r\ -\ ) (5) 

The same equivalence can be shown for the variable A = 
m\Silm,2 — Si /mi). That is, since 

q = (l-5m)/(l + 6m), (6) 

we find 

A = {S 2 ~S 1 )-2.Sm.(S 2 +Si)+2.Sm 2 .(S 2 -S 1 ) + ..., (7) 



and hence A can be reexpressed in terms of the alter- 
native spin variables. However, as we saw in |25j . the 
dominant contribution to the recoil is actually described 
most simply in terms of the variables r\ and A (in partic- 
ular, by the expression ?y 2 AjJ. 

Assuming a Taylor expansion of a given order of ex- 
pansion (o) one can count how many terms (n) there 
should be if we have (v) independent variables using the 
following multinomial formula [28] 



(o + v- 1)! 
o\(v- 1)! ' 



(8) 



However, only certain combinations of variables are al- 
lowed. In order to take into account the correct combina- 
tions of variables for each component of the recoil velocity 
at a given order, we consider the symmetry properties 
summarized in Tabic ||J The possible terms to a given 
expansion order in spin (i.e., products of S and A) are 
summarized in Tables HlllIVI Note that in these terms co- 
efficients can depend on higher, even and odd, powers of 
8m respectively. The terms in Tables |III| and |IV| are all 
multiplied by fitting coefficients. 



Interestingly, for the (5m—independent terms, we can 
obtain the spin-dependence of the in-plane recoil from 
the spin-dependence of the out-of-plane recoil via 

V± = Vf|[A_L <— > A,,], (9) 

while for the 5m— dependent terms that are odd powers 
in the spin variables , we have 

V x {Sm) = Sm.V\\ [S ± <— > S {] ]. (10) 

On the other hand, for terms proportional to even powers 
of the spin variables there are extra terms not present in 
Vj|. 

In addition, functionally, the terms proportional to 5m 
in Vji can be obtained from the 5m- independent terms in 

V\\ by 

V\\(5m) = 5m.V\\[S± «— > AJ, (11) 

while for odd powers of the spins only, the i5m-dependent 
terms in V± can be obtained from the 5m-independcnt 
terms via 

VI (5m) = 6m.V±[S\\ < — > Ay]. (12) 

No such correspondence holds for the for the 5m- 
dependent terms with odd powers in the spin for V±. 

B. The equal-mass case 

Using the above properties we find 16 terms up to 
fourth-order in the spin that contribute to the off-plane 



TABLE II: Number of possible terms at a given order of 
expansion (with respect to 5* or A. Here 1 indicates terms 
present even in the equal-mass limit (and/or proportional to 
even powers of 8m) and 5m indicates terms proportional to 
5m to odd powers. 



Order 


0th 


0th 


1st 


1st 


2nd 


2nd 


3rd 


3rd 


4th 


4th 


total 


mass 


1 


5m 


1 


5m 


1 


8m 


1 


5m 


1 


8m 


All 


Vj| 








1 


1 


2 


2 


5 


5 


8 


8 


32 


Vx 





1 


1 


1 


2 


4 


5 


5 


8 


11 


38 


Total 





1 


2 


2 


4 


6 


10 


10 


16 


19 


70 



recoil velocity 



Vj| = A± + A ± .S\\ + A t S ± + 
A ± .A| + Ai + A ± .5f + 
A ± .Si + A|| .S±.S\\ + 
Ai.Sj, +A||.Si + A!.Sj_ + A_ L .£f + 

A||.s , j_.s| + A||.a3_.0 l + 

A ± .5||.A| + A ± .5||.5i + --- (13) 
We can regroup all these terms (assuming we can col- 
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TABLE III: Parameter dependence at each order of expansion for the off-plane recoil. Here 1 indicates terms present even in 
the equal-mass limit (and/or proportional to even powers of Sm) and Sm indicates terms proportional to 8m to odd powers. 



II 


Oth order 


1 





Sm 





Vj, 


1st order 


1 


A ± 


Sm 


Sx 


Vj| 


2nd order 


1 


A x .5|| +A|,.5x 


Sm 


Ax-Ay + S±.S\\ 


V\\ 


3rd order 


1 


A,, .Sx.Sy + A ± .Sf + Ax.Ajj + Ai + A±. S'i 


Sm 


S x .A| + Sx.Sf +A x .A||.5||+5x.Ai+5i 


Vj, 


4th order 


1 


Sx.Ajj + A x .Sjf + Ax.S||.Ajj + 5x-A||.5jf + Al.Sy + SI .Ay + A ± .S X .A|| + Ax.S x .S|| 


Sm 


Ax-Ay + S ± .SI + 5x.5y.A| + Ax-Ay.Sf 4- A x .Ay 4- S x .Sy + A X .S X .S|| 4- S X A X .A|| 



TABLE IV: Parameter dependence at each order of expansion for the in-plane recoil. Here 1 indicates terms present even in 
the equal-mass limit (and/or proportional to even powers of Sm) and Sm indicates terms proportional to Sm to odd powers. 





Oth order 


1 





Sm 


1 


v ± 


1st order 


1 


Ay 


Sm 




Vx 


2nd order 


1 


Ay. Sy 4- Ax.Sx 


Sm 


Ajj + sf + Al + S'i 


Vx 


3rd order 


1 


Ax.Sx.Sn 4- Ay.Sjf 4- Ay.A x + Ajj + Ay.S x 4- Ajj 


Sm 


S,| .Ajj + S,|. S'i + A ± . S X .A|| +S||.Ai + Sf 


v ± 


4th order 


1 


Sx.Ai + A||.Sjf 4- A||.Sy.AI 4- Sx.A x .Sff 4- Af.Sy 4- S X .A X 4- Ajj.S x .A ± + A||.SI.S|| 


Sm 


Ax.A||.Sx.Sy 4- Ai + A| + Si 4- S[[ 4- A X .A| + Al.Sl 4- Al.Sf 4- Ajj.S x 4- Ajj.Sjf 4- Sl-Sjf 



etc, should really be expressed as CiA-fi4-C2A-TOS 2 4-- ■ • , 
i.e., not only are there fitting constants Ci, C%, • • • , but 
each coefficient also has its own angular dependence. We 
will return to this issue in Sec. EH 



III. NUMERICAL RELATIVITY TECHNIQUES 

For the black-hole binary (BHB) data presented here, 
both BHs have the same mass, but have different spins. 
We use the TwoPUNCTURES thorn [53] to generate ini- 
tial puncture data [30] for the BHB simulations described 
below. These data are characterized by mass parameters 
m pi/2i momenta Pi/%, spins S1/2, and coordinate loca- 
tions X1/2 of each hole. We obtain parameters for the 
location, momentum, and spin of each BH using the 2.5 
PN quasicircular parameters. Here we choose to nor- 
malize the PN initial data such that the total Arnowitt- 



lect all J_ terms) in the following form 

V l{ = A x .(l4-A 2 ± 4----)-(l + A|4----)- 

.(1 + Si + • • • ).(1 + S|| + Sf + Sf + • • • ) + 

+5x.A||.(l + Ai + •■■).(! + Af + •••)• 

.(l + Sl4----)-( 1 + S||+Sf + -..). (14) 

The 10 terms directly proportional to the cos(cp) are 
those linear to the subindex J_ 

^1°^= A_i_.(l + Ajj + ..-).(l + S|| + 4-- •■) + 
+ Sx.Ay.(H- A| 4- •••)•(! + S\\ +Sf ] +---) (15) 

This (symbolic) expression is the one we will use in this 
paper. 

There is a subtlety in the above expansion. Because S± 
and A x are vector quantities, terms like A±(l + S z + - ■ ■ ), 
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Deser-Misner (ADM) energy is 1M. We obtain param- 
eters rripi/2 using an iterative procedure in order to ob- 
tain a system where the two BHs have the same mass 
and the total ADM energy is 1M. This iterative proce- 
dure is most efficient when the horizon masses and ADM 
energy can be obtained from the initial data alone. For 
highly-spinning BHs (a = S/m 2 > 0.9), a relatively large 
amount of energy lies outside the BH. This energy is even- 
tually absorbed, changing the mass of the BH substan- 
tially (see, e.g. [31]). We therefore limit the spin of the 
BHs to a =< 0.8. 

We evolve these BHB data sets using the LazEv [32] 
implementation of the moving puncture approach [2] [3] 
with the conformal function W = ^Jx = cxp(— 20) sug- 
gested by Ref. [33J. For the runs presented here, we use 
centered, eighth-order finite differencing in space |34j and 
a fourth-order Runge Kutta time integrator. (Note that 
we do not upwind the advection terms.) 

Our code uses the EinsteinToolkit [35] [3J] / Cac- 
tus [37] / Carpet [35] infrastructure. The Carpet 
mesh refinement driver to provide a "moving boxes" style 
of mesh refinement. In this approach refined grids of 
fixed size are arranged about the coordinate centers of 
both holes. The Carpet code then moves these fine 
grids about the computational domain by following the 
trajectories of the two BHs. 

We obtain accurate, convergent waveforms and horizon 
parameters by evolving this system in conjunction with a 
modified 1+log lapse and a modified Gamma-driver shift 
condition [2] |3SJ HO], and an initial lapse a(t = 0) = 
2/(l + ipg L ), where ipBL is the Brill-Lindquist conformal 
factor and is given by 



BL 



= 1 



X/(2|F-r1|), 



where ri is the coordinate location of puncture i. The 
lapse and shift are evolved with 

(d t - li l di)a = -2aK, (16a) 
d t (3 a = (3/4)f Q - 7#\ (16b) 



momentum, and angular momentum, in terms of the ra- 
diative Weyl Scalar ip^, using the formulas provided in 
Refs. [?3J HI]. However, rather than using the full ipi, 
we decompose it into £ and m modes and solve for the 
radiated linear momentum, dropping terms with t > 5. 
The formulas in Refs. j43j|44] are valid at r = oo. We ex- 
tract the radiated energy-momentum at finite radius and 
extrapolate to r = oo using both linear and quadratic 
extrapolations. We use the difference of these two ex- 
trapolations as a measure of the error. 

Both the remnant parameter variation, and the varia- 
tion in the extrapolation to infinity of the radiation un- 
derestimate the actual errors in the quantity of interest. 
However, because quantities like the total radiated en- 
ergy can be obtained from either extrapolations of ipi or, 
quite independently, from the remnant BHs mass, the 
difference between these two is a reasonable estimate for 
the actual error. 

Our empirical formula will depend on the spins mea- 
sured with respect to the orbital plane at merger. In 
Ref [2 5) we described a procedure for determining an ap- 
proximate plane. This is based on locating three fiducial 
points on the BHBs trajectory r+, rg, and r_, where r+ 
is the point where r(t) (r(t) is the orbital separation) 
reaches its maximum, r_ is the point where r(t) reaches 
its minimum, and r*o is the point between the two where 
r(t) = 0. These three points can then be used to define 
an approximate merger plane (see Fig. [2]). We then need 
to rotate each trajectory such that the infall directions 
all align (as much as possible). This is accomplished by 
rotating the system, keeping the merger plane's orienta- 
tion fixed, such that the vector r + — r is aligned with 
the y axis. The azimuthal angle ip, described below, is 
measured in this rotated frame. 

Our motivation for defining the orbital plane "at 
merger" is the observation that most of the recoil is gen- 
erated near (and slightly after) merger. For example, 
Fig. [3] shows the recoil imparted to the remnant BH for 
the N45PH30 configuration. As seen in the plot, all but 
16% of the recoil is generated "post-merger." 



where we use r\ = 2 for all simulations presented below. 

We use AHFinderDirect [41] to locate apparent 
horizons. We measure the magnitude of the horizon 
spin using the isolated horizon (IH) algorithm detailed 
in Ref. [32]. Note that once we have the horizon spin, 
we can calculate the horizon mass via the Christodoulou 
formula 

™H = y/m* r + S%/(4m* I ) i (17) 

where m; rr = WAJ (I6n) and A is the surface area of 
the horizon, and Sh is the spin angular momentum of 
the BH (in units of M 2 ). In the tables below, we use 
the variation in the measured horizon irreducible mass 
and spin during the simulation as a measure of the error 
in these quantities. We measure radiated energy, linear 



IV. SIMULATIONS 

In this paper we consider four families of equal-mass, 
precessing, BHB configurations, which we will denote by 
S, K, L, and N. Initial data parameters are given in Ta- 
bles [x] and [Xl]( found in Appendix All. These configura- 



tions are characterized by the spins of the two BHs on the 
initial slice. While we have chosen different orientations 
of the spins, the magnitudes of the intrinsic spin of each 
BH in the binary ai — 0.8, with the exception of the N 
configurations, where the the first BH has spin a\ = 0.8 
and the second is nonspinning. For the S configurations, 
Si = —S2. While for the K configurations S\ z = —S2Z 
but S\ x — £>2x an d Siy — S^y. The L configurations 
have the spin of BH1 entirely in the orbital plane, while 
the spin of BH2 is perpendicular to the plane, and fi- 
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FIG. 3: A plot of the absolute values of P z (t) (the momentum 
imparted to the remnant) and dP z /dt(t) for NTH45PH30 con- 
figuration. The vertical line represents the approximate time 
of merger. 




FIG. 2: Finding the orbital plane near merger. The upper plot 
shows the orbital separation r(t) versus time. The inset shows 
r(t) near merger and r(t) (rescaled by 200 for clarity). The 
points f + , r*o, and r_ correspond to the times where r is max- 
imized, zero, and minimized (denoted with arrows here). The 
plot below shows the trajectory, the points f+,fo,f- (large 
red dots) and the "merger" plane. 



S2=0 





FIG. 4: The N configuration. These configurations differ from 
the hangup kick configurations in that BH2 is nonspinning. 
Numerical evolutions preserve the N configurations only ap- 
proximately. 



nally the N configurations have BH1 spinning and BH2 
nonspinning. We use the notation zTHxxxPHyyy, where 
z is N, S, K, or L, xxx gives the inclination angle 9 of 
spin of BH1 and yyy gives the orientation of the spin of 
BH1 in the initial orbital plane, i.e. the azimuthal angle 
4>. In order to fit the resulting recoils, we found that we 
needed at least six azimuthal configurations in the inter- 
val [0, 180°) for each 9 configurations. This is due to the 
fact that we need to separate contributions due to cos 30 
from contributions due to cos <f> (by symmetry, the recoil 
out of the plane cannot contain terms of the form cos n<f) 
if n is even). 

In all cases, the computational domain extended to 
±400M, with a coarsest resolution of h = AM at the 
outer boundary. We used 9 levels of refinement, centered 
on each puncture, with radii 200, 100, 50, 20, 10, 5, 2, 



0.6, respectively. 

We chose these configurations for two main reasons, 
first each family of configurations can be described by a 
single azimuthal angle parameter <j> and a single polar an- 
gle 9, and they activate different terms in our ansatz for 
the recoil. The former is necessary in order to reduce the 
computational costs. In general, four angular parame- 
ters are required in order to describe the spins at merger 
(two polar and two azimuthal). In order to model the 
polar and azimuthal dependence, we would need at least 
6x6x6x6 = 1296 simulations (per choice of spin magni- 
tude, per choice of mass ratio). By reducing the dimen- 
sionality to two, we only need 36 simulations (per family, 
per a, per q). This reduction only works, however, if the 
two parameter family of initial data, maps in a straight- 
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FIG. 5: The S configuration. These configuration differ from 
the hangup kick configuration in that Si z = — S2 Z (and hence 
Si = —02), initially. Numerical evolutions preserve the S 
configurations only approximately. 
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FIG. 6: The K configuration. These can be thought of as 
a modification of the S configurations. There Si z = — S2 Z , 
while Sin: = S2X and Si y = S2 y , initially. 
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TABLE V: S, K, and L configuration angles. Here (j>i is the 
angle between the in-plane component of Si for configuration 
PHO and the corresponding PHXX configuration, while 02 is 
defined using S2, 4>a using A, and 4>s using S. These angles 
agree to within (3 — 15)°, which justifies our using a single 
angle in our fitting formula for this configuration. For the S 
configurations, 4>s is ill-defined since Sj_ is very small. 



Conf 


<t>i 


02 


4>a 


d>s 


STH45PH0 














STH45PH30 


33.3518 


30.9202 


32.0924 




STH45PH60 


61.7152 


58.2499 


59.9753 




STH45PH90 


91.4444 


88.5921 


90.0456 




STH45PH120 


116.271 


115.641 


115.997 




STH45PH90 


143.981 


144.881 


144.46 




KTH45PH0 














KTH45PH30 


28.6011 


31.0559 


31.4676 


28.8129 


KTH45PH60 


39.3732 


41.6249 


42.4836 


39.2248 


KTH45PH90 


58.7044 


59.2109 


61.8581 


56.8317 


KTH45PH105 


75.1457 


72.9151 


76.7906 


71.6609 


KTH45PH120 


95.7513 


90.7761 


95.0918 


91.3502 


KTH45PH135 


114.209 


108.253 


112.123 


110.072 


KTH45PH150 


135.705 


130.366 


133.203 


132.59 


KTH45PH165 


156.52 


153.385 


154.838 


154.882 


KTH22.5PH0 














KTH22.5PH30 


29.9 


23.9 


26.8 


25.12 


KTH22.5PH60 


52.1 


42.1 


45.7 


46.3 


KTH22.5PH90 


71.4 


59.8 


62.8 


71.4 


KTH22.5PH120 


101.2 


92.8 


92.9 


102.2 


KTH22.5PH150 


142.5 


143.4 


139.9 


147.5 


LPHO 














LPH30 


17.9 


20.8 


17.0 


24.7 


LPH60 


35.0 


36.3 


32.0 


45.3 


LPH90 


55.2 


52.3 


49.8 


65.3 


LPH120 


87.6 


80.0 


81.4 


93.0 


LPH150 


139.9 


130.8 


136.1 


138.4 



forward way to a 2 parameter family of configurations 
at merger. In particular, we need the final configuration 
to be describable by a single azimuthal and polar angle. 
We note that this is not the case in general, and that we 
used PN simulations to tests the stability of various con- 
figurations. For the N configuration, the mapping to a 
single azimuthal angle is automatic because only one BH 
is spinning, for the other configurations, we verify that 
the configuration can be described by a single angle by 
comparing four different measurements of the azimuthal 
angle. The results are displayed in Table [Vj which shows 
that, to within about 3-15 degrees, the configurations are 
describable by a single angle. 



RESULTS 



FIG. 7: The L configuration. The L configuration is a mod- 
ification of the N configuration, where Si is aligned with the 
orbital angular momentum z and rather than having S2 = 0, 
S2 is varied initially in the orbital plane. 



Results from these 83 simulations are given in the ta- 
bles in Appendix |A2| below. In Tables |XII| and |XII 
we give the remnant BH mass and spin, as measured 
using the IH formalism, while in Tables |XIV| and |XV| 
we give the radiated energy, angular momentum, and 



recoil, as calculated from the waveform extracted at 
60M, 70M, ■■■ , 100M and then extrapolated to infin- 
ity. We compare these two independent measures in Ta- 
bles |XVT| and prvTl] Finally, in Tables |XVIII| ancfXTX) 
we give the spin of each BH near merger, and the final 
remnant recoil, in a rotated frame aligned with averaged 
orbital angular momentum at merger. For completeness, 
we also show in Table |XX| the value of Aj_ , S± , A z , and 
S z corresponding to the BH spin in Table XVIH 



In order to analyze the results of the present simula- 
tions, we use the techniques developed in [35]. Briefly, 
we rotate each configuration such that the trajectories 
near merger overlap. We then calculate the spins in this 
rotated frame. This is done separately for each family of 
constant 9 per configuration type (S,L,K,N). The angle 
<p is then defined to be the angle (at merger) between 
spin of BH1 (the BH originally located on the positive 
x axis) for a given PHyyy configuration and the spin of 
BH1 in the corresponding PHO configuration. Note that, 
for a given family of fixed spin and spin inclination angle 
9, the angle ip and <p differ by a constant, which can be 
absorbed in the fitting constants fa and fa. We then fit 
the recoil in these sub-families to the form 



Vi cos(<^ - fa) ± V 3 cos(3<y9 - 3fa). 



(18) 



Our tests indicate that V\ can be obtained accurately 
with 6 choices of the initial fa angles. For example, a fit 
all the NTH45PHyyy gi ves V Y = 1349.0±9.7 if we include 
all 12 angles (see Table |XVIlij ), and V x = 1346 ±22 if we 
include 6 angles. 

In all cases, V3 is much smaller than V\. Results from 
these fits are given in Table |VI| and Fig. [8j We note 
that there are additional approximations inherent in this 
procedure. To demonstrate this, consider the formula 



V g = A ± (A + BS Z 



■) + S ± A z (D + ES z 



), (19) 



where A, B,C, D, E are fitting constants. Even when 
considering "symmetric" configurations like S, K, L, and 
N, where where S and A cannot rotate independently, 
each term in Eq. ( 19 ) may be maximized at different az- 



imuthal angles, and the formula should really be written 
as 



V g = AA ■ h + BA ■ hiS z 
+DS ■ h 2 A z + ES ■ h 3 A z S z - 



(20) 



where hi are unit vectors in the orbital plane. If we 
make the additional assumption that the coefficient A 



dominates this expression, then Eq. (20) can be approx- 
imated by 



V z = A-h (A + Bcos(tf 1 )S z ) 
h A z (D cos(# 2 ) + Ecos($ 3 )S z ), 



(21) 



where $j is the angle between hi and ho. There will be 
terms proportional to sin$i, but they will be 0(1/ A), 
which we will assume to be small enough to ignore. 



TABLE VI: A fit of the recoil for each family of PHYYY 
configurations to the form V\\ = Vi cos(<y9 — fa) + V3 cos(3<^ — 
3^3). Note how the K configurations, which started with 
Ax = 0, evolved to configurations with large Ai. 



CONF 



NTH15 

NTH30 

NTH45 

NTH60 

NTH120 

NTH135 

NTH165 



STH45 



KTH45 
KTH22.5 



Vi 



Vz 



fa 



539.34 ±2.5 
1002 ± 12 
1349.0 ±9.7 
1542 ± 11 
1199 ± 13 
927.5 ±6.4 
312.9 ±6.4 



33.2 ±2.3 
43 ± 13 
52 ± 12 
34 ± 11 
37 ± 12 
35.6 ±6.7 
11.6 ±6.2 



141.96 ±0.24 
126.42 ±0.71 
82.50 ±0.58 
20.83 ±0.47 
292.79 ± 0.54 
226.90 ± 0.43 
213.4 ± 1.2 



297.1 ± 1.3 
260.3 ±5.3 
337.0 ±3.9 

269.2 ±6.8 
139.6 ±5.9 

311.3 ±3.6 
189 ± 11 



2020 ± 19 50 ± 19 291.40 ± 0.56 342.3 ± 7.3 



2227 ± 12 
1731 ± 25 



195 ± 12 
130 ± 23 



217.33 ±0.32 
164.49 ±0.75 



155.0 ± 1.3 
100.9 ±3.4 



3014 ±21 145 ± 18 331.30 ±0.36 263.50 ±2.9 




3D 100 




60 • BO. ■ 100 120 140 



20 40 B60 Sb 100 120 [40 



FIG. 8: A fit of the recoil for the NTH45PHyyy (left) and 
NTH135PHyyy (right) configurations to V = Vx cos(<£ — </>i) 
(blue-dotted) and V = Vicoa{ip — fa)- s t-vzcos(^(p — ^fa) (red- 
solid), as well as the residuals (blue circles and red squares, 
respectively). The remaining scatter is due to the fact that 
Si_ and S z vary from configuration to configuration within 
the NTH45PHyyy and NTH135PHyyy families. 



If, in addition, we assume that the angles i)i do not 
vary significantly between different configurations, then 
we can replace the fitting constants in Eq. ( plj ) by the 
product of the constant with the corresponding cos^. 



We can then interpret V\ from Eq. ( 18 ) as the V\ = 
\A ± \(A+BS Z )+\S ± \\A Z \(D+ES Z ), where B — Bcostfi, 
etc. Ultimately, we justify all our approximations by test- 
ing the resulting formula using several different families 
of configurations. 

The N configurations are the only ones with families of 
different 9 (apart from the two angles for the K config- 
uration). For these families, we fit the data in Tabic VI 
to the form 

V\— C\Ot sin 9 ± C 2 a 2 sin 9 cos 9 ± 

C 3 a 3 sin 9 cos 2 9 ± VLmg, (22) 

where Phang is t ne hangup kick [171 118) . which has the 



9 



100 ISO 




III 120 140 i 



1500 



1000 



0.5 



1.0 



2.0 



2.5 



3.0 



FIG. 9: Fitting of the KTH22.5PHyyy (left) and 
KTH45PHyyy (right) configurations to V = Vi cos(<£ — </>i) 
(blue-dotted) and V = V\ cos(ip — (f>i ) + v 3 cos(3<p — 3^3) (red- 
solid), as well as the residuals (blue circles and red squares, 
respectively). Note how strong the higher-order contributions 
are compared to the other configurations. 



10 



TABLE VII: Fit of the NTH data (see Table |VIf to the form 
Vi = Cia sin 9 + C 2 a 2 sin 9 cos 9 + C 3 a 3 sin 8 cos 2 9 + Viang- 
Corrections to the hangup kick are significant at the (assumed) 
a 2 and higher orders. 



Coeff. 


Correction to hangup kick 


hangup kick term 


Ci 


41.0 ±46.8 


3677.76 


c 2 


1220 ± 96 


2481.21 


c 3 


2104 ±669 


1792.45 


c 4 





1506.52 



form 

V hang = 3677.76a sin 9 + 2481. 21a 2 sin 9 cos 9 + 
1792.45a 3 sin6<cos 2 6>± 
1506.52a 4 sinfl cos 3 9. (23) 

For consistency with our conventions in j!7l 118] , we define 
a here to be a x /2. We note that S\\/M 2 = (l/2)a and 
Ay/M 2 = a. We also note that, because a < 1/2 here, 
the extrapolation to a = 1 is more severe than in the 
original hangup kick configuration. Results from these 
fits are shown in Table IVTT1 



VI. MODELING THE RECOIL VELOCITY 

When considering a generic formula for the recoil, we 
need to take into account discrete symmetry operations, 
such as parity (see Ref. [371 US)- Motivated by PN, we 
construct a formula for the recoil using combinations of 
S = (Si + §2), A = m(S l 2 /TO2 — Si/ m i)j an d Sm — 
mi — TO2, where Sifi is measured during the late inspiral, 
and m = mi + mi- The formula must be symmetric 
under exchange of labels lf»2. Of particular interest to 
us will be the quantities Aj^ = A • h/m 2 , where n is a 



0.5 



2.0 



2.5 



FIG. 10: A fit of the NTHXXX families, with residuals, and 
comparison with the hangup kick formula (black dotted) and 
supcrkick (more symmetrical, blue-dotted curve). An excess 
over both formulas at small angles is apparent, while a slight 
deficit with respect to the hangup kick is apparent at large 
angles. 



unit vector in the orbital plane (at merger), S± = S ■ h, 
and the z components of S and A. 

Considering the equal-mass case, and up to cubic order 
in the spin, the most general formula for the out-of-plane 
recoil is 



v; 



A ± (A + BS Z + CS 2 Z ) + 
S ± A Z (D + ES Z ) + 
FAj_Sj_ + 
GAj_A 2 + 
HAl, 



(24) 



where A — H are fitting constants. The first line is part 
of the hangup kick recoil, the second and third are new 
contributions. The term proportional to G is small (if G 
were big, then the S configuration recoils would be signif- 
icantly different from the hangup kick prediction) . The H 
term is small, as we saw by evolving superkick configura- 
tion in |20) . Motivated by the hangup kick results, where 
the series A + BS Z + CS Z had similarly larger values of 
A, B, and C, we will assume at this point that the term 
proportional to E is larger than the term proportional to 
F. 



We therefore interpret the additional terms in Eq. ( 22 ) 
as being proportional to powers of A z and S z and 
corrections to the hangup kick formula have the form 
2C2S_\_A z +AC^S_i_A z S z . Then, if we assume the same rf 
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TABLE VIII: Comparison of Vi as fit from the current data and the predictions of the superkick, hangup kick, and cross kick 
(new) formulas. Note, there is an ambiguity in the sign of the cross kick correction for the S configuration (see text). 



CONF 


S ± /M 2 


A ± /M 2 


S z /M 2 


A z /M 2 


Vi 


Sup. 


Hang. 


Cross. 


NTH 15 


0.046 ±0.004 


0.092 ±0.008 


0.196 ±0.001 


-0.392 ±0.002 


539.4 ±2.3 


339.746 


463.256 


537 


NTH30 


0.090 ±0.007 


0.179 ±0.013 


0.179 ±0.003 


-0.358 ±0.007 


1002 ± 12 


658.497 


871.282 


998 


NTH45 


0.126 ±0.008 


0.252 ±0.015 


0.155 ±0.006 


-0.311 ±0.013 


1349.0 ±9.7 


926.499 


1176.76 


1322 


IN 1 ribU 


0.161 ± 0.0073 


0.323 ± 0.015 


0.118 ± 0.010 


r\ not; i r\ non 

—0.235 ± 0.020 


1542 ± 11 


1186.7 


1413.2 


1543 


NTH120 


0.184 ± 0.004 


0.368 ± 0.008 


-0.077 ±0.010 


0.154 ±0.021 


1199 ± 13 


1355.8 


1279 


1178 


NTH135 


0.154 ±0.006 


0.308 ±0.011 


-0.128 ±0.007 


0.256 ±0.013 


927.5 ±6.4 


1134.46 


967.015 


913 


NTH165 


0.059 ±0.003 


0.118 ±0.007 


-0.193 ±0.002 


0.386 ± 0.004 


312.9 ± 6.4 


434.141 


342.312 


324 


KTH45 


0.276 ± 0.002 


0.497 ± .028 


0.054 ±0.021 


-0.277 ±0.048 


2227 ± 12 


1826 


1970 


2191 


KTH22.5 


0.149 ± .003 


0.400 ± 0.037 


0.021 ±0.008 


-0.626 ±0.025 


1731 ± 25 


1470 


1512 


1756 


L 


0.173 ±0.016 


0.551 ±.006 


0.227 ±0.013 


0.103 ±0.051 


3014 ± 21 


2026 


2928 


3006 


STH45 


0.011 ± 0.004 


0.552 ± 0.004 


0.005 ± 0.003 


-0.5760 ±0.0015 


2020 ± 19 


2030.15 


2044.94 


2060* 



mass ratio dependence, our ansatz for the z component 
of the generic recoil becomes 



Vz 



= Aj_ (3677.76 + 2 x 2481.21S* 

1-4 x 1792.45S'2 ± 8 x 1506.525^) 
-S±A Z (2C 2 +AC 3 S Z ). (25) 



We refer to this new contribution to the recoil, which has 
the form S±A Z (2C2 ± iC^Sz), as the cross kick (since 
S±A 2 can be expressed as z ■ S x (h x A), where n is a 
unit vector in the xy plane). The coefficients Ci and C3 
were determined using only the N configurations. 

We then verify Eq. ( 25 ) , in the equal-mass limit, by 



comparing the predictions from the new formula with the 
maximum recoil obtained from the S, K, and L configura- 



tions. Our results are given in Table VIII| The table com- 
pares the measured value of V\ for each family with the 
predictions of the superkick, hangup kick, and cross kick. 
In all cases, except S, the cross kick provides the most ac- 
curate prediction for V\ . The results from the K configu- 
rations arc particularly interesting since the measured re- 
coils and the cross kick predictions are both ~ 200 km s _1 
larger than the hangup kick prediction (a 10-16% effect). 
For the S configurations, there is an ambiguity in the 
sign of the cross kick kick. This is due to the fact that 
the cross kick correction lies in the same direction as the 
hangup kick. Here, however, S z = 0.005 ± 0.003 and the 
small hangup kick correction may have the wrong sign. 
If we assume S z is really zero or slightly negative, we get 
find that the cross kick prediction is the most accurate 
(with a prediction of 2015 km s -1 ). 



VII. DISCUSSION 

The discovery that the hangup effect contributes sig- 
nificantly to the gravitational recoil of merging black hole 
binaries |17j implies that nonlinear spin couplings are 
crucial in describing those recoils. Nonlinear couplings 
come in a variety of combinations, as described in Sec. [Til 
In order to evaluate which of those terms produce the 
largest contributions to the total recoil, we performed a 
large set of new simulations. These 83 simulations of pre- 
cessing BHBs allowed us to confirm the relevance of the 
hangup kick effect in more generic runs, discover another 
similarly important term that we named cross kick that 
appears in precessing binaries, and gives more accurate 
predictions for other families of BHB configurations. 

While not as dramatic as the hangup kick effect, the 
cross kick may prove to be very important in the non- 
equal-mass regime. To help elucidate how this new con- 
tribution affects the recoil (for a given mass ratio), we 
plot the maximum recoil for configuration with a given 
mass ratio and with both BHs maximally spinning. As 



shown in Fig. 11 the cross kick enhances the recoil (up 
to 600 km s _1 ) in the moderate mass-ratio range. 

To see how the cross kick contribution to the recoil af- 
fects the net probabilities for large recoils, we revisit the 
case of the supermassive BH binary with spins aligned 
via hot and cold accretion [18]. Briefly, we consider a 
set of 10 million binaries chosen randomly with a spin- 
magnitude distribution and spin inclination angle distri- 
bution taken from [18] . and a mass ratio distribution 
taken from 46-48J. We assume a uniform distribution 
of spin directions in the equatorial plane. We find an 
increased probability of large recoils (V > 2000 km s _1 ) 
by a factor of ~ 2 (see Table IX However, to generate 
these pr oba bilities we used the assumption that all terms 
in Eq. (25 1 scale with the mass ratio q as I677 2 , where 
i] = q/(l + q) 2 . This is a strong assumption that we will 
revisit in an upcoming paper. Additionally, we did not 
take into account new nonlinear terms of the form 8m ■ ■ ■ 
that may also prove to be important. 

As an aside, we note that the distribution of azimuthal 
orientations near merger may appear to be nonuniform. 
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r\ (symmetric mass ratio) 

FIG. 11: The maximum recoil velocity predicted by the 
hangup kick and cross kick formulas for BHBs with a given 
mass ratio and maximal spin. The inset shows the difference 
between these two velocities versus symmetric mass ratio. 



TABLE IX: Comparison between the predicted probabilities 
for a recoil in a given range as from the hangup kick and cross 
kick formulas for hot (top) and cold (bottom) accretion. 



Range 


P (cross) 


P (cross obs) 


P(hang) 


P(hang obs) 


0-500 


77.000% 


91.301% 


80.871% 


93.210% 


500-100 


15.564% 


6.903% 


13.843% 


5.623% 


1000-2000 


6.930% 


1.741% 


5.046% 


1.143% 


2000-3000 


0.498% 


0.055% 


0.237% 


0.025% 


3000-4000 


0.007% 


3.5 • 10 _4 % 


0.003% 


10" 4 % 


0-500 


91.193% 


97.765% 


93.657% 


98.522% 


500-100 


7.974% 


2.114% 


5.919% 


1.423% 


1000-2000 


0.832% 


0.120% 


0.423% 


0.055% 


2000-3000 


0.002% 


1.3- 10~ 4 


4.7- 10~ 4 % 


% 


3000-4000 


0% 


0% 


0% 


0% 



This is actually most pronounced for the a = 0.9 in 
our original hangup kick paper |18j (see Fig. 5 there). 
However, this skew appears to be actually due to vary- 
ing eccentricity, which leads to different inspiral times 
for different starting azimuthal configurations. To help 
confirm that there is not, in fact, a strong preference 
for any particular azimuthal angles, we evolved as set 
of 360 superkick configurations (with the spins aligned 
along 4> = 0°, 1°, • • • , 359°), using 3.5 PN, from a sepa- 
ration of 10M to 3M (note, we are not concerned with 
the accuracy of PN at 3M, rather, if there is any signif- 
icant effect predicted by PN). The distribution of final 
azimuthal configurations was flat, with no strong pref- 
erence or clumping. A plot of final versus the initial 
azimuthal angle <f> is shown in Fig. [12] There is a small 
sinusoidal effect at the level of 4 parts in 1000. 

By examining the N configuration (which have non- 
zero values for Sn, S±, Am, and Aj_), we found a new 
nonlinear term that amplifies the recoil. We verified this 




0.994 - 
0.992 - 



FIG. 12: A plot of d0finai/d</>init versus time for a set of 360 
binaries in a superkick configuration. The initial separation 
is ~ 10M, while the final separation is ~ 3M. The effect is 
only 4 parts in 1000. 

new effect be examining several other configurations (S, 
K, L). Since the N configurations are generic, in that all 
relevant spin parameters are non-trivial, it appears to be 
the case that there is no other large nonlinear contribu- 
tion to the recoil for equal-mass BHBs (such a contri- 
bution would have to vanish when one the BH is non- 
spinning). On the other hand, the unequal-mass regime, 
which is the subject of major research effort by the au- 
thors, promises to hold many new surprises. 
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Appendix A: Progenitor and Remnant Parameters 

The tables in this appendix provide useful information 
for modeling remnant properties and how they relate to 
the configuration of the progenitor BHB. These results 
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can be used to improve, both quantitatively and qualita- 
tively, the empirical formulas that describe the remnant 
mass and spin [2H US]. They can also be used in the 
construction of alternative formulas to model recoil ve- 
locities. Although we have been able to accurately model 



the off-plane recoil and make predictions that fit new 
runs, the formulas become increasingly complex to model 
higher nonlinear terms in the spins and one can seek a 
simpler, more compact formulation of the remnant recoil. 



1. Initial Data parameters 
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TABLE X: Initial data parameters. In all cases the puncture masses were chosen such that the total ADM mass of the binary 
was 1.0± 10 _6 M. Here the punctures are located at (£1,2, 0, 0) with momenta ±(0,p, 0) and spins Si = (S x , S y , S z )- For the N 
configurations S2 = 0. The approximate initial eccentricities a, eccentricities measured over the last orbit e/, and the number 
of orbits N, are also given. 



CONF 


rripi/M 




xi/M x 2 /M 


p/M 


S x /M 2 


Sy/M 2 


S z /M 2 


mm 


m H 2 


NTH15PH0 





O A77C 

307753 





a 77 A a 

487740 


4 


077900 


—4 


307566 





1 A70AO 

107293 


000000 


053084 


0.198112 





r A700 1 

507321 





C A A 7 C A 

504754 


NTH15PH30 





A 77/"* A 

307762 





A 77 A 1 

487741 


4 


077900 


-4 


307566 





1 A7A A O 

107293 


—0.026542 


0.045972 


0.198112 





P" A7000 

507322 





504756 


NTH15PH60 





307779 





487745 


4 


077900 


4 


307566 





1 A70AO 

107293 


—0.045972 


0.026542 


0.198112 





r A700 1 

507321 





C A A 7 E A 

504759 


NTH15PH90 





A77C") 

307786 





A O 77 A 7 

487747 


4 


A^'TAAA 

077900 


4 


307566 





1 A7A A O 

107293 


—0.053084 


0.000000 


0.198112 





f A70 A 1 

507321 





504761 


NTH15PH120 





307777 





A O 77 A E 

487745 


4 


077900 


-4 


307566 





1 A70 A O 

107293 


—0.045972 


—0.026542 


0.198112 





r A7Q 1 A 

507319 





E A /I 7 E A 

504759 


NTH15PH150 





307763 





A O 77 A 1 

487741 


4 


077900 


-4 


307566 





1 A70AO 

107293 


—0.026542 


—0.045972 


0.198112 





r A7QO 1 

507321 





504755 


NTH30PH0 





307652 





a ^7^7 a r* 

487705 


4 


098641 


-4 


304988 





1 A A A 1 

107491 


000000 


102538 


177601 





507291 





r* A J ^7 A A 

504702 


NTH30PH30 

1 1 1 -1-lOU-L llOU 





307686 





A O ^7^7 1 A 

487710 


4 


098641 


-4 


304988 





1 A A A 1 

107491 


— 051269 


088800 


1 77601 





P" A ^7A A O 

507293 





r* A A ^7 A 

504708 


NTH30PH60 





a 77 a a 

307744 





A O r7ITA P" 

487725 


4 


098641 


-4 


OA A A O O 

304988 





1 A A A 1 

107491 


— 088800 


051269 


177601 





P" A ^7 A A 1 

507291 





r* a j ^7 a 1 

504721 


NTH30PH90 





A 7777^7 7 

307775 





J ^7^70 1 

487731 


4 


098641 


-4 


OA J A O O 

304988 





1 A A A 1 

107491 


—0 1 02538 


000000 


1 77601 

U.l 1 1 UU 1 





P" A ^7 A O O 

507288 





r" A A ^7 A ^7 

504727 


NTH30PH120 





/~i P717 ^ 1- 

307745 





,1 O 77A f 

487725 


4 


AAO C A -\ 

098641 


-4 


OA A AO O 

304988 





1 A7 A A 1 

107491 


— 088800 


—0 051269 


177601 

l_/«ll 1 UU 1 





P* A 7A O 7 

507287 





f A ,1 rr A A 

504720 


NTH30PH150 





307685 





487711 


4 


098641 


-4 


304988 





107491 


— 051269 


—0 088800 


1 77601 

U.l 1 1 UU 1 





507289 





504707 







307512 





487654 


4 


131521 


-4 


300556 





107804 


000000 


144983 


144983 





507241 





504625 


NTH45PH10 





307522 





487655 


4 


131521 


-4 


300556 





107804 


—0 0251 76 


142780 


144983 





507244 





504626 


NTH45PH20 





307544 





487659 


4 


131521 


-4 


300556 





107804 


—0 04Q587 


1 36239 


144983 





507244 





504631 


NTH45PH30 





307584 





487664 


4 


131521 


-4 


300556 





107804 


—0 072491 


1 25^59 


144983 





507246 





504635 


NTH45PH40 





307615 





A 7/ 1 7/" 1 

487676 


4 


131521 


-4 


O A A f f /' 

300556 





1 A70 A A 

107804 


— 093193 


1 1 1063 


144983 

U. 1 1 A CO u 





P* A 7A A O 

507242 





f r\ a a 7 

504647 


NTH45PH45 





307636 





J A 7/"* A 1 

487681 


4 


131521 


-4 


O A A 7 t"~ J ' 

300556 





1 A fro A /I 

107804 


— 102518 


102518 


144983 

Ui 1 1 1 JUU 





P" A ^7 A A A 

507242 





504651 


NTH45PH50 





307658 





,1 O 7/" 1 O P* 

487685 


4 


131521 


-4 


O A A f f 

300556 





1 A70 A A 

107804 


— 111063 


093193 


144983 

Ui 1 1 1 - 7 V 1 * ) 





P* A 7A A 1 

507241 





504655 


NTH45PH60 





307700 





487693 


4 


131521 


-4 


O A A 7 f"~ J ' 

300556 





1 A fro A /I 

107804 


— 125559 


072491 


144983 

U ■ 1 1 1 JUU 





507240 





504663 


NTH45PH75 

IN 1 11 1 U 





307746 





487703 


4 


131521 


-4 


300556 





107804 


— 140042 


037524 


144983 

U. 1 1 1 JUu 





507238 





504671 


NTH45PH90 

1\ 1 ll^Ul J. J. U 





r\77i~ 7 

307757 





A 77ao 

487708 


4 


131521 


-4 


O A A f f /' 

300556 





1 A70 A /I 

107804 


— 144983 


000000 


144983 

u. 1 njuu 





P* A 7A O O 

507233 





r* a < A7r 

504675 


NTH45PH975 





O A77E C 

307755 





/I O 77A7 

487707 


4 


131521 


-4 


300556 





1 A70 A /I 

107804 


—0.143742 


—0 018924 


144983 





r a 70 O O 

507233 





E A A a 70 

504673 


NTH45PH120 





307699 





487694 


4 


131521 


-4 


oaac rrn 

300556 





1 A70 A /I 

107804 


-0.125559 


-0.072491 


0.144983 





r A 70 OA 

507230 





E A A Ct O A 

504660 


NTH45PH1125 





0A77O0 

307723 





/I O 77AH 

487700 


4 


131521 


-4 


OAAr cr C 

300556 





1 A70 A /I 

107804 


-0.133947 


—0.055482 


0.144983 





p» A700A 

507229 





e a /i a 00 

504666 


NTH45PH1275 





307668 





/i 7r* 

487688 


4 


131521 


-4 


a a rr c c 

300556 





1 A70 A /I 

107804 


-0.115022 


-0.088260 


0.144983 





rr A700A 

507229 





504655 


NTH45PH150 





O A7T O 1 

307581 





487665 


4 


131521 


-4 


OAAf" f /' 

300556 





1 A70 A A 

107804 


—0.072491 


—0.125559 


0.144983 





P* A 7A O 

507236 





504634 


NTH60PH0 





a 70 7 i 

307374 





487595 


4 


i7/ii 7 a 

174170 


-4 


OA A 1 i~\f 

294196 





108211 


000000 


1 77524 


1 02493 

U. 1U^t:i7U 





P* A71 OA 

507180 





504534 


NTH60PH30 





n7 a a c 

307465 





/i 7r* 1 c 

487616 


4 


17/11 7fl 

174170 


-4 


OA A 1 C\G 





108211 


—0 088762 


153740 

V7»_LOO 1 IV' 


102493 





CT A71 OA 

507184 





504556 


NTH60PH60 





307656 





A 7r* e z? 

487656 


4 


i7/ii 7n 

174170 


-4 


OA A 1 C\G 





108211 


-0.153740 


0.088762 


0.102493 





CT A71 OA 

507180 





E A A E AO 

504592 


NTH60PH90 





n77 /i 

307743 





487679 


4 


i7/ii 7n 

174170 


-4 


OA A 1 C\G 





108211 


—0.177524 


0.000000 


0.102493 





CT A71 7 A 

507170 





504611 


NTH60PH120 





307652 





/I O *7C E 7 

487657 


4 


17/11 7fl 

174170 


-4 


OA A 1 A£? 





108211 


-0.153740 


-0.088762 


0.102493 





CT A71 O A 

507164 





CT c\ A con 

504589 


NTH60PH150 





O fl7 /I It 

307465 





/I O *7d 1 7 

487617 


A 

4 


17/11 7fl 

174170 


4 


OA A 1 A£? 

294iyD 





108211 


—0.088762 


—0.153740 


0.102493 





CT A71 d A 

507169 





504552 


NTH120PH0 





a 7r\7 

307974 





A O O O OA 

488280 


4 


1 7 

661137 


4 


P* A A A 1 A 

542010 





1 A I - A O A 

105032 


000000 


177129 


—0.102265 





506644 





f f\ A A 1 /I 

504014 


NTH120PH30 





308059 





,1 O O O A 

488303 


4 


1 7 

661137 


-4 


P* /I A A 1 A 

542010 





1 A I - A O A 

105032 


— 088564 


153398 


—0.102265 





506646 





r* a ^ a 7 

504037 


NTH120PH60 










/I O O O A O 


4 


etc 1 1 7 

661137 


-4 


r A OA 1 A 

54^U1(J 





1 A cr A OO 

105032 


—0.153398 


0.088564 


—0.102265 





pr C\CC A O 





E A A A70 

504073 


NTH120PH90 





no 007 

308327 





488364 


4 


661137 


-4 


r /I OA 1 A 

54^U1(J 





1 A f. A O O 

105032 


—0.177129 


0.000000 


—0.102265 





rACCO /I 

506634 





E A A AAA 

504090 


NTH120PH120 





noo a 1 





A O O O A O 


4 


ca 1 1 07 

661137 


-4 


r /I OA 1 A 

54^U1(J 





1 A ti A O O 

105032 


-0.153398 


-0.088564 


-0.102265 





506630 





e a a nr 

504068 


NTH120PH150 





308058 





,1 /"i 
488303 


4 


1 7 

661137 


-4 


P* A OA 1 A 

542010 





1 A I - A O O 

105032 


—0.088564 


—0.153398 


—0.102265 





506635 





r* A 1 A O A 

504032 


NTH135PH0 

li 1 111UU1 11U 





308101 





488288 


4 


701321 


-4 


532307 





105392 


000000 


144599 


—0 1 44599 





506603 





504004 


NTH1 ^PFRO 





308159 





488302 


4 


701321 


-4 


532307 





105392 


— D799QQ 









506602 





504019 


NTH135PH60 





308281 





488329 


4 


701321 


-4 


532307 





105392 


-0.125226 


0.072299 


-0.144599 





506600 





504043 


NTH135PH90 





308342 





488343 


4 


701321 


-4 


532307 





105392 


-0.144599 


0.000000 


-0.144599 





506597 





504053 


NTH135PH120 





308282 





488329 


4 


701321 


-4 


532307 





105392 


-0.125226 


-0.072299 


-0.144599 





506591 





504039 


NTH135PH150 





308161 





488301 


4 


701321 


-4 


532307 





105392 


-0.072299 


-0.125226 


-0.144599 





506595 





504014 


NTH165PH0 





308524 





488582 


4 


854936 


-4 


625367 





104087 


0.000000 


0.052895 


-0.197408 





506462 





503917 


NTH165PH30 





308529 





488584 


4 


854936 


-4 


625367 





104087 


-0.026448 


0.045809 


-0.197408 





506461 





503920 


NTH165PH60 





308549 





488587 


4 


854936 


-4 


625367 





104087 


-0.045809 


0.026448 


-0.197408 





506463 





503922 


NTH165PH90 





308553 





488590 


4 


854936 


-4 


625367 





104087 


-0.052895 


0.000000 


-0.197408 





506461 





503925 


NTH165PH120 





308545 





488588 


4 


854936 


-4 


625367 





104087 


-0.045809 


-0.026448 


-0.197408 





506461 





503923 


NTH165PH150 





308529 





488584 


4 


854936 


-4 


625367 





104087 


-0.026448 


-0.045809 


-0.197408 





506460 





503919 



4.5 



4 -0().()()3 



4.5 



^0.02 

4 0.003 



3.5 



o c;0.02 



q r0.02 
J-OO.OOS 
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TABLE XI: Initial data parameters. In all cases the puncture masses were chosen such that the total ADM mass of the binary 
was 1.0 ± 10~ 6 M. Here the punctures are located at (xi,2,0,0) with momenta ±(0,p,0) and spins Si = (S x , S y , S z ). For the 
S configurations the second BH spin is given by S2 = —Si, while for the K configurations it is given by S2 = {S x , S y , —S z ). 
Finally, for the L configurations, S2 = (0, 0, \Si\). The approximate initial eccentricities e», eccentricities measured over the 
last orbit e/, and the number of orbits N, are also given. 



CONF 



m p i/M m p2 /M xi/M x 2 /M p/M S x /M 2 S y /M 2 S z /M 2 m H i m H 2 



STH45PH0 





302923 





303004 


4 


194252 


-4 


534000 





107508 


0.000000 


0.144813 





144813 





505422 





505411 


^•5o!oo4 


STH45PH30 





303051 





303137 


4 


194252 


-4 


534000 





107508 


-0 


072406 


0.125412 





144813 





505446 





505436 




STH45PH60 





303313 





303398 


4 


194252 


-4 


534000 





107508 


-0 


125412 


0.072406 





144813 





505482 





505474 




STH45PH90 





303442 





303531 


4 


194252 


-4 


534000 





107508 


-0 


144813 


0.000000 





144813 





505496 





505488 




STH45PH120 





303312 





303400 


4 


194252 


-4 


534000 





107508 


-0 


125412 


-0.072406 





144813 





505468 





505461 




STH45PH150 





303051 





303137 


4 


194252 


-4 


534000 





107508 


-0 


072406 


-0.125412 





144813 





505433 





505424 




KTH45PH0 





303010 





303092 


4 


187859 


-4 


527100 





107491 


0.000000 


0.144822 





144822 





505460 





505454 


«J-°0.004 


KTH45PH30 





303090 





303173 


4 


187859 


-4 


527100 





107491 


-0 


072411 


0.125419 





144822 





505469 





505459 




KTH45PH60 





303250 





303335 


4 


187859 


-4 


527100 





107491 


-0 


125419 


0.072411 





144822 





505474 





505464 




KTH45PH90 





303330 





303416 


4 


187859 


-4 


527100 





107491 


-0 


144822 


0.000000 





144822 





505471 





505465 




KTH45PH105 





303310 





303393 


4 


187859 


-4 


527100 





107491 


-0 


139887 


-0.037483 





144822 





505466 





505462 




KTH45PH120 





303249 





303336 


4 


187859 


-4 


527100 





107491 


-0 


125419 


-0.072411 





144822 





505460 





505458 




KTH45PH135 





303171 





303254 


4 


187859 


-4 


527100 





107491 


-0 


102404 


-0.102404 





144822 





505457 





505454 




KTH45PH150 





303091 





303173 


4 


187859 


-4 


527100 





107491 


-0 


072411 


-0.125419 





144822 





505456 





505452 




KTH45PH165 





303028 





303118 


4 


187859 


-4 


527100 





107491 


-0 


037483 


-0.139887 





144822 





505456 





505454 




KTH45PH30 





303090 





303173 


4 


187859 


-4 


527100 





107491 


-0 


072411 


0.125419 





144822 





505469 





505459 




KTH45PH60 





303250 





303335 


4 


187859 


-4 


527100 





107491 


-0 


125419 


0.072411 





144822 





505474 





505464 




KTH45PH90 





303330 





303416 


4 


187859 


-4 


527100 





107491 


-0 


144822 


0.000000 





144822 





505471 





505465 




KTH45PH120 





303249 





303336 


4 


187859 


-4 


527100 





107491 


-0 


125419 


-0.072411 





144822 





505460 





505458 




KTH45PH150 





303091 





303173 


4 


187859 


-4 


527100 





107491 


-0 


072411 


-0.125419 





144822 





505456 





505452 





KTH22.5PH0 

KTH22.5PH30 

KTH22.5PH60 

KTH22.5PH90 

KTH22.5PH120 

KTH22.5PH150 



0.303316 
0.303338 
0.303385 
0.303406 
0.303384 
0.303336 



0.303422 
0.303447 
0.303494 
0.303520 
0.303495 
0.303449 



4.133612 
4.133612 
4.133612 
4.133612 
4.133612 
4.133612 



-4.576616 
-4.576616 
-4.576616 
-4.576616 
-4.576616 
-4.576616 



0.107486 
0.107486 
0.107486 
0.107486 
0.107486 
0.107486 



0.000000 

-0.039189 

-0.067878 

-0.078379 

-0.067878 

-0.039189 



0.078379 

0.067878 

0.039189 

0.000000 

-0.039189 

-0.067878 



0.189223 
0.189223 
0.189223 
0.189223 
0.189223 
0.189223 



0.505509 
0.505512 
0.505514 
0.505512 
0.505510 
0.505508 



0.505497 
0.505500 
0.505502 
0.505503 
0.505501 
0.505499 



0.505627 5lHos 

0.505652 

0.505703 

0.505727 

0.505699 

0.505649 



LPH0 

LPH30 

LPH60 

LPH90 

LPH120 

LPH150 



0.302983 
0.303110 
0.303359 
0.303478 
0.303361 
0.303111 



0.303118 
0.303193 
0.303348 
0.303430 
0.303345 
0.303193 



4.400426 
4.400426 
4.400426 
4.400426 
4.400426 
4.400426 



-4.164556 
-4.164556 
-4.164556 
-4.164556 
-4.164556 
-4.164556 



0.105611 
0.105611 
0.105611 
0.105611 
0.105611 
0.105611 



0.000000 

-0.102508 

-0.177549 

-0.205015 

-0.177549 

-0.102508 



0.205015 

0.177549 

0.102508 

0.000000 

-0.102508 

-0.177549 



0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 



0.505728 
0.505737 
0.505732 
0.505722 
0.505720 
0.505722 
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2. Remnant Properties and Radiated Energy-Momentum 



TABLE XII: Remnant horizon properties using the IH formalism. Quoted errors are calculated from the variation of IH 
quantities with time. See Tables [XVI| and |XVII| for more realistic measures of the true error. 



CONF 


m H /M 












a 














S x /M 








Sy/M 






S z /M 










NTH15PH0 





949002 


± 


2 


5 


10" 


b 





800400 


± 


2 


9 


10" 


b 


-0.0059188 ± 


5 


10" 


a 


0.0229940 ±5 


• IO - " 







72045 


± 


3 


10" 


4 


NTH15PH30 





949096 


± 


2 


5 


10" 


-5 





800619 


± 


3 


1 


10" 


5 


-0.0087750 ± 


5 


10" 


7 


0.0327727 ± 2 


• 10" 6 







72039 


± 


4 


10" 


-5 


NTH15PH60 





949177 


± 


2 


5 


10" 


-5 





800721 


± 


2 


8 


10" 


5 


-0.0193710 ± 


1 


10" 


6 


0.0200562 ± 1 


• 10" 6 







72086 


± 


4 


10" 


-5 


NTH15PH90 





949169 


± 


2 


6 


10" 


-5 





800619 


± 


3 


4 


10" 


5 


-0.0247011 ± 


5 


10" 


4 


0.0037424 ± 2 


• 10" 3 







72086 


± 


4 


10" 


-5 


NTH15PH120 





949076 


± 


2 


6 


10" 


-5 





800397 


± 


3 


3 


10" 


5 


-0.0231286 ± 


1 


10" 


2 


-0.0123739 ± 


7- 10" 


3 





72048 


± 


4 


10" 


-4 


NTH15PH150 





948988 


± 


2 


5 


10" 


-5 





800289 


± 


2 


9 


10" 


5 


-0.0271884 ± 


1 


10" 


2 


-0.0324823 ± 


7- 10" 


3 





71948 


± 


1 


10" 


-3 


NTH30PH0 





950352 


± 


1 


1 


10" 


-5 





793089 


± 


8 


7 


10" 


6 


0.0102968 ± 3 


• io- 7 




0.0606032 ± 2 


• 10~ 6 







71365 


± 


2 


10" 


-5 


NTH30PH30 





950660 


± 


1 





10" 


-5 





793625 


± 


1 


1 


10" 


5 


-0.0433712 ± 


6 


10" 


7 


0.0612195 ±8 


■IO" 7 







71331 


± 


1 


10" 


-5 


NTH30PH60 





950722 


± 


1 


1 


10" 


-5 





793366 


± 


1 


9 


10" 


5 


-0.0469123 ± 


7 


10" 


7 


0.0326284 ±5 


■IO" 7 







71482 


± 


2 


10" 


-5 


NTH30PH90 





950488 


± 


1 


3 


10" 


-5 





792563 


± 


1 


2 


10" 


5 


-0.0599863 ± 


1 


10" 


6 


0.0046845 ± 2 


■IO" 7 







71349 


± 


2 


10" 


-5 


NTH30PH120 





950142 


± 


1 


2 


10" 


-5 





791909 


± 


9 


2 


10" 


6 


-0.0680924 ± 


2 


10" 


6 


-0.0250607 ± 


8 • 10" 


7 





71122 


± 


3 


10" 


-5 


NTH30PH150 





950063 


± 


9 


9 


10" 


-6 





792167 


± 


1 


2 


10" 


5 


-0.0279134 ± 


1 


10" 


6 


-0.0566108 ± 


2 • 10" 


6 





71223 


± 


3 


10" 


-5 


NTH45PH0 





952965 


± 


4 


2 


10" 


-6 





781452 


± 


1 


8 


10" 


5 


0.0052411 ±2 


•HT 7 




0.0862119 ±2 


• 10" 6 







70439 


± 


2 


10" 


-5 


NTH45PH10 





952963 


± 


5 





10" 


-6 





781286 


± 


2 


4 


10" 


5 


-0.0068446 ± 


3 


10" 


7 


0.1028876 ±3 


• IO" 6 







70198 


± 


3 


10" 


-5 


NTH45PH20 





952898 


± 


5 





10" 


-6 





780915 


± 


2 


7 


10" 


5 


-0.0280650 ± 


9 


10" 


7 


0.0982404 ± 3 


• 10" 6 







70168 


± 


3 


10" 


-5 


NTH45PH30 





952788 


± 


3 


9 


10" 


-6 





780408 


± 


2 


5 


10" 


5 


-0.0396175 ± 


1 


10" 


6 


0.0727714 ±2 


• IO" 6 







70360 


± 


3 


10" 


-5 


NTH45PH40 





952657 


± 


3 


7 


10" 


-6 





779872 


± 


2 


4 


10" 


5 


-0.0573442 ± 


2 


10" 


6 


0.0647666 ± 2 


• 10" 6 







70247 


± 


2 


10" 


-5 


NTH45PH45 





952594 


± 


4 





10" 


-6 





779622 


± 


2 


8 


10" 


5 


-0.0611883 ± 


2 


10" 


6 


0.0605063 ± 2 


• IO" 6 







70220 


± 


3 


10" 


-5 


NTH45PH50 





952517 


± 


5 


1 


10" 


-6 





779337 


± 


3 


3 


10" 


5 


-0.0622867 ± 


2 


10" 


6 


0.0533615 ±2 


• KT 6 







70231 


± 


3 


10" 


-5 


NTH45PH60 





952371 


± 


5 


5 


10" 


-6 





778824 


± 


3 


5 


10" 


5 


-0.0819302 ± 


4 


10" 


6 


0.0635694 ±3 


• IO" 6 







69875 


± 


4 


10" 


-5 


NTH45PH75 





952129 


± 


5 


7 


10" 


-6 





778112 


± 


4 


3 


10" 


5 


-0.0837902 ± 


4 


10" 


6 


0.0188236 ±9 


■IO" 7 







70015 


± 


4 


10" 


-5 


NTH45PH90 





951927 


± 


8 


6 


10" 


-6 





777677 


± 


5 





10" 


5 


-0.0926429 ± 


5 


10" 


6 


0.0144582 ± 9 


■IO" 7 







69844 


± 


5 


10" 


-5 


NTH45PH975 





951866 


± 


6 


8 


10" 


-6 





777638 


± 


5 


1 


10" 


5 


-0.0940831 ± 


6 


10" 


6 


-0.0000789 ± 


5 • 10" 


8 





69827 


± 


5 


10" 


-5 


NTH45PH120 





951917 


± 


7 


5 


10" 


-6 





778302 


± 


5 


8 


10" 


5 


-0.0863067 ± 


6 


10" 


6 


-0.0311843 ± 


2 ■ 10" 


6 





69926 


± 


5 


10" 


-5 


NTH45PH1125 





951859 


± 


6 





10" 


-6 





777949 


± 


5 


2 


10" 


5 


-0.0884714 ± 


5 


10" 


(i 


-0.0238097 ± 


1 • 10" 


6 





69887 


± 


5 


10" 


-5 


NTH45PH1275 





952012 


± 


5 


3 


10" 


-6 





778733 


± 


4 


8 


10" 


5 


-0.0813575 ± 


5 


10" 


6 


-0.0424656 ± 


3 • 10" 


6 





69979 


± 


5 


10" 


-5 


NTH45PH150 





952488 


± 


4 


1 


10" 


-6 





780369 


± 


2 


3 


10" 


5 


-0.0509898 ± 


2 


10" 


6 


-0.0707933 ± 


2 • 10" 


6 





70258 


± 


3 


10" 


-5 


NTH60PH0 





954548 


± 


9 


4 


10" 


-6 





759903 


± 


3 


6 


10" 


5 


0.0120181 ±6 


■IO' 7 




0.1059528 ±5 


• IO" 6 







68413 


± 


3 


10" 


-5 


NTH60PH30 





954117 


± 


8 


1 


10" 


-6 





758701 


± 


3 


8 


10" 


5 


-0.0608939 ± 


2 


10" 


(i 


0.1016092 ±4 


• 10~ 6 







68044 


± 


3 


10" 


-5 


NTH60PH60 





954263 


± 


8 


3 


10" 


-6 





758808 


± 


4 





10" 


5 


-0.1015464 ± 


4 


10" 


6 


0.0580071 ± 2 


• IO" 6 







68102 


± 


3 


10" 


-5 


NTH60PH90 





954825 


± 


9 


1 


10" 


-6 





760276 


± 


1 


9 


10" 


5 


-0.1254426 ± 


4 


10" 


6 


-0.0000259 ± 


5 • 10" 


8 





68169 


± 


3 


10" 


-5 


NTH60PH120 





955661 


± 


5 


2 


10" 


-6 





762668 


± 


2 


3 


10" 


5 


-0.0923975 ± 


2 


10" 


6 


-0.0521792 ± 


1 ■ 10" 


6 





68841 


± 


2 


10" 


-5 


NTH60PH150 





955595 


± 


5 


8 


10" 


-6 





762609 


± 


1 


9 


10" 


5 


-0.0525913 ± 


9 


10" 


7 


-0.1010084 ± 


2 • 10" 


6 





68701 


± 


2 


10" 


-5 


NTH120PH0 





965664 


± 


3 


6 


10" 


-6 





644178 


± 


8 


9 


10" 


6 


-0.0196764 ± 


3 


10" 


7 


0.1049745 ±2 


• IO" 6 







59113 


± 


1 


10" 


-5 


NTH120PH30 





965986 


± 


3 


5 


10" 


-6 





644558 


± 


8 


9 


10" 


6 


-0.0503997 ± 


7 


10" 


7 


0.0898719 ± 1 


• 10~ 6 







59256 


± 


1 


10" 


-5 


NTH120PH60 





965901 


± 


3 


7 


10" 


-6 





643317 


± 


8 


3 


10" 


6 


-0.1022275 ± 


1 


10" 


6 


0.0432199 ± 6 


■IO" 7 







58984 


± 


1 


10" 


5 


NTH120PH90 





965596 


± 


4 


6 


10" 


-6 





641770 


± 


8 


5 


10" 


6 


-0.1096112 ± 


1 


10" 


(i 


-0.0101321 ± 


2 • 10" 


7 





58816 


± 


1 


10" 


-5 


NTH120PH120 





965027 


± 


4 


3 


10" 


-6 





640323 


± 


1 


2 


10" 


5 


-0.0745858 ± 


1 


10" 


6 


-0.0501520 ± 


8 ■ 10" 


7 





58951 


± 


1 


10" 


-5 


NTH120PH150 


o 


965083 


± 


3 


5 


10" 


-6 





641861 


± 


1 


1 


10" 


5 


-0.0473661 ± 


8 


10" 


7 


-0.0969227 ± 


2 • 10" 


6 


o 


58801 


± 


1 


10" 


-5 


NTH135PH0 





967523 


± 


4 


8 


10" 


-6 





611243 


± 


1 


4 


10" 


5 


-0.0111685 ± 


2 


10" 


7 


0.0886077 ± 2 


• IO" 6 







56517 


± 


2 


10" 


-5 


NTH135PH30 





967277 


± 


3 


8 


10" 


-6 





609924 


± 


7 


9 


10" 


6 


-0.0536160 ± 


8 


10" 


7 


0.0687996 ± 1 


• 10" 6 







56395 


± 


1 


10" 


-5 


NTH135PH60 





967256 


± 


4 





10" 


-6 





609695 


± 


9 


3 


10" 


6 


-0.0842716 ± 


1 


10" 


6 


0.0381523 ± 6 


•IO" 7 







56287 


± 


1 


10" 


-5 


NTH135PH90 





967402 


± 


4 


4 


10" 


-6 





610218 


± 


1 





10" 


5 


-0.0831299 ± 


1 


10" 


6 


-0.0206110 ± 


4- 10" 


7 





56462 


± 


1 


10" 


-5 


NTH135PH120 





967739 


± 


3 


5 


10" 


-6 





611715 


± 


8 


3 


10" 


6 


-0.0587466 ± 


9 


10" 


7 


-0.0630778 ± 


1 ■ 10" 


6 





56636 


± 


1 


10" 


-5 


NTH135PH150 





967820 


± 


3 


7 


10" 


-6 





612348 


± 


1 





10" 


5 


-0.0352962 ± 


6 


10" 


7 


-0.0865760 ± 


1 • 10" 


6 





56590 


± 


1 


10" 


-5 


NTH165PH0 





969939 


± 


4 


2 


10" 


-6 





566054 


± 


1 





10" 


5 


-0.0143038 ± 


3 


10" 


7 


0.0238197 ±4 


•IO" 7 







53181 


± 


1 


10" 


-5 


NTH165PH30 





969916 


± 


3 


5 


10" 


-6 





565949 


± 


9 


5 


10" 


6 


-0.0157146 ± 


3 


10" 


7 


0.0275530 ± 5 


•IO" 7 







53146 


± 


1 


10" 


-5 


NTH165PH60 





969934 


± 


3 


9 


10" 


-6 





565965 


± 


9 


7 


10" 


6 


-0.0247061 ± 


4 


10" 


7 


0.0148809 ± 3 


•IO" 7 







53166 


± 


1 


10" 


-5 


NTH165PH90 





969973 


± 


4 


5 


10" 


-6 





566113 


± 


9 


6 


10" 


6 


-0.0399334 ± 


6 


10" 


7 


-0.0092692 ± 


2 • 10" 


7 





53105 


± 


1 


10" 


-5 


NTH165PH120 





969995 


± 


4 


3 


10" 


-6 





566233 


± 


8 


5 


10" 


6 


-0.0282134 ± 


5 


10" 


7 


-0.0310082 ± 


5- 10 - 


7 





53111 


± 


1 


10" 


-5 


NTH165PH150 





969978 


± 


5 


1 


10" 


-6 





566194 


± 


8 


8 


10" 


6 


-0.0058152 ± 


1 


10" 


7 


-0.0245145 ± 


4- 10" 


7 





53211 


± 


1 


10" 


-5 
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TABLE XIII: Remnant horizon properties using the IH formalism (continued). Quoted errors are calculated from the variation 
of IH quantities with time. See Tables [XVI| and |XVII| for more realistic measures of the true error. 



CONF 


m H /M 






a 






S x /M 






Sy/M 




S z /M 






STH45PH0 


0.960386 ±1.2 


10" 




0.681328 ±5.5 


10" 




-0.0029258 ±2 


10" 




-0.0112825 ±7-10" 




0.62831 ± 4 


10" 




STH45PH30 


0.961370 ± 1.1 


10" 




0.683685 ± 3.8 


10" 




-0.0049041 ± 2 


10" 




-0.0106043 ±4- 10" 




0.63178 ± 3 


10" 




STH45PH60 


0.961447 ± 1.2 


10" 




0.682989 ±4.3 


10" 




-0.0059768 ± 3 


10" 




-0.0059344 ±3-10" 




0.63128 ± 3 


10" 




STH45PH90 


0.960884 ± 1.3 


10" 




0.680984 ± 5.8 


10" 




0.0043894 ± 3 • 10" 7 




0.0118551 ±7 - 10" 7 




0.62862 ± 4 


10" 




STH45PH120 


0.960038 ± 1.5 


10" 




0.678902 ± 5.9 


10" 




-0.0020680 ± 2 


10" 




-0.0083257±6- 10" 




0.62567 ±4 


10" 




STH45PH150 


0.959700 ± 1.6 


10" 




0.678717 ±6.1 


10" 




-0.0006975 ± 9 


10" 




0.0130117±8- 10" 7 




0.62498 ± 4 


10" 




KTH45PH0 


0.959306 ± 4.6 


10" 




0.733285 ±1.5 


10" 




0.0096842±3- 10"' 




0.1670485±4 - 10" b 




0.65374 ± 2 


10" 


b 


KTH45PH30 


0.959031 ± 4.6 


10" 


g 


0.731242 ± 7.0 


10" 


-6 


-0.0824495 ± 1 


10" 


Q 


0.1443256 ± 2 • 10~ 6 




0.65169 ± 1 


10" 


-5 


KTH45PH60 


0.958963 ± 4.3 


10" 


g 


0.729211 ± 3.8 


10" 




-0.1581502 ±7 


10" 


Q 


0.0892345 ± 4 • 10" 6 




0.64554 ± 3 


10" 




KTH45PH90 


0.958292 ± 6.6 


10" 


g 


0.727374 ±1.3 


10" 




-0.1681021 ±4 


10" 


g 


-0.0020425 ±1-10" 


j 


0.64646 ± 2 


10" 




KTH45PH105 


0.957553 ± 1.3 


10" 




0.725788 ±4.7 


10" 




-0.1589683 ±6 


10" 


g 


-0.0343141 ±1-10" 




0.64530 ± 3 


10" 




KTH45PH120 


0.956969 ± 7.6 


10" 


g 


0.724630 ± 2.1 


10" 




-0.1604927 ±3 


10" 


g 


-0.0888209 ±2-10" 


g 


0.63775 ± 2 


10" 




KTH45PH135 


0.957068 ± 3.3 


10" 


g 


0.725770 ± 1.3 


10" 




-0.1359932 ±3 


10" 


g 


-0.1290310 ±3-10" 




0.63781 ± 2 


10" 




KTH45PH150 


0.958008 ± 3.9 


10" 


g 


0.729194 ± 3.7 


10" 




-0.0820269 ±4 


10" 


g 


-0.1513764 ± 7 ■ 10" 


g 


0.64671 ± 3 


10" 


-5 


KTH45PH165 


0.958990 ±4.4 


10" 


g 


0.732376 ± 2.2 


10" 




-0.0389777 ± 1 


10" 


g 


-0.1725190 ±5-10" 


(j 


0.64990 ± 2 


10" 


_5 


KTH22.5PH0 


0.959881 ±5.2 


10" 




0.695101 ±2.4 


10" 




-0.0103440 ± 3 


10" 




0.0872087±2- 10" b 




0.63440 ± 2 


10" 




KTH22.5PH30 


0.959651 ±6.5 


10" 




0.694469 ± 2.8 


10" 




-0.0551628 ±2 


10" 




0.0763874 ± 2 • 10" 6 




0.63258 ± 2 


10" 




KTH22.5PH60 


0.959845 ± 5.0 


10" 




0.694970 ± 1.8 


10" 




-0.0860895 ± 2 


10" 




0.0396870 ± 8 • 10" 7 




0.63322 ± 2 


10" 




KTH22.5PH90 


0.960412 ±3.7 


10" 




0.696387 ± 1.6 


10" 




-0.0946340 ± 2 


10" 




-0.0085426 ±2-10" 




0.63527 ± 2 


10" 








1 n~ 


6 




1 n~ 


-5 


-0.0768168 ±2 


10" 


6 


n n' : !799i 4- q . i n~ 


7 


n fi^S^Q 4- 9 
u.uoojy in z 


1 D" 


-5 


KTH22.5PH150 


0.960575 ± 4.8 


10" 


6 


0.696820 ± 1.4 


10" 


-5 


-0.0395611 ±6 


10" 


7 


-0.0746323 ±1-10" 


6 


0.63739 ± 1 


10" 


-5 


LPHO 


0.943756 ± 2.4 


10" 


5 


0.815716 ±2.8 


10" 


4 


-0.0000661 ± 8 


10" 


8 


0.1262529±5- 10" b 




0.71548 ±3 


10" 


4 


LPH30 


0.943158 ±9.2 


10" 


6 


0.814545 ±4.7 


10" 


-5 


-0.0688881 ±5 


10" 


6 


0.1112364±8- 10" 6 




0.71266 ±6 


10" 


-5 


LPH60 


0.943272 ±2.1 


10" 


5 


0.814734 ±3.1 


10" 


-4 


-0.1118922 ±5 


10" 


5 


0.0791998±3- 10" 5 




0.71184 ± 3 


10" 


-4 


LPH90 


0.944022 ± 2.2 


10" 


5 


0.816420 ±3.4 


10" 


-4 


-0.1400266 ±6 


10" 


5 


-0.0082355 ±4- 10" 


6 


0.71393 ±3 


10" 


-4 


LPH120 


0.945434 ± 2.6 


10" 


5 


0.819694 ±4.2 


10" 


-4 


-0.1281539 ±7 


10" 


5 


-0.0585698 ±3-10" 


5 


0.71900 ± 4 


10" 


-4 


LPH150 


0.945209 ± 1.8 


10" 


5 


0.818582 ± 1.7 


10" 


-4 


-0.0745095 ± 2 


10" 


5 


-0.1087893±3- 10" 


5 


0.71935 ± 2 


10" 


-4 
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TABLE XIV: Radiated mass, angular momentum, and the remnant recoil (in original frame) as calculated from 1^4. Errors 
quoted are from differences between to extrapolation to r = 00. See Tables [XVI| and |XVII| for more accurate measurement of 
the error. 



CONF 


<5M rad 








SJ X 








5Jy 




8.L 








V x 


Vy 


V, 


NTH15PH0 





0494 


± 





0001 


-0.0015 


± 





0008 


0.0199 ±0 


0023 





3690 


± 





0023 


75 ±22 


113 ±5 


-456 ± 2 


NTH15PH30 





0493 


± 





0001 


-0.0090 


± 





0004 


0.0154 ±0 


0005 





3687 


± 





0024 


75 ±21 


121 ±5 


-248 ± 2 


NTH15PH60 





0492 


± 





0001 


-0.0140 


± 





0016 


0.0068 ± 


0015 





3685 


± 





0025 


67 ± 21 


126 ±5 


21 ± 1 


NTH15PH90 





0492 


± 





0001 


-0.0153 


± 





0022 


-0.0035 ± 


0.0029 





3690 


± 





0021 


60 ±21 


125 ±4 


286 ± 1 


NTH15PH120 





0493 


± 





0001 


-0.0126 


± 





0022 


-0.0129 ± 


0.0037 





3689 


± 





0025 


59 ±21 


117±5 


482 ±2 


NTH15PH150 





0494 


± 





0001 


-0.0065 


± 





0017 


-0.0189 ± 


0.0035 





3692 


± 





0023 


68 ±22 


111 ± 5 


546 ± 2 


NTH30PH0 





0481 


± 





0001 


-0.0040 


± 





0021 


0.0377 ±0 


0036 





3588 


± 





0028 


129 ± 21 


41 ±6 


-548 ± 3 


NTH30PH30 





0478 


± 





0001 


-0.0194 


± 





0009 


0.0285 ± 


0003 





3576 


± 





0032 


116 ± 19 


71 ±5 


-58 ±2 


NTH30PH60 





0477 


± 





0001 


-0.0288 


± 





0015 


0.0120 ±0 


0039 





3575 


± 





0035 


85 ± 16 


81 ±6 


417 ± 1 


NTH30PH90 





0480 


± 





0001 


-0.0324 


± 





0013 


-0.0075 ± 


0.0060 





3587 


± 





0036 


64 ± 16 


61 ±7 


809 ±2 


NTH30PH120 





0483 


± 





0001 


-0.0271 


± 





0010 


-0.0252 ± 


0.0070 





3602 


± 





0031 


74 ± 18 


30 ±8 


1019 ±3 


NTH30PH150 





0483 


± 





0001 


-0.0140 


± 





0011 


-0.0362 ± 


0.0059 





3601 


± 





0029 


108 ± 20 


19 ±8 


927 ±3 


NTH45PH0 





0456 


± 





0001 


-0.0075 


± 





0041 


0.0507 ±0 


0015 





3410 


± 





0039 


103 ± 12 


-6 ±4 


190 ± 1 


NTH45PH10 





0456 


± 





0001 


-0.0167 


± 





0051 


0.0481 ± 


0000 





3409 


± 





0041 


85 ±11 


10 ±3 


425 ± 1 


NTH45PH20 





0456 


± 





0001 


-0.0258 


± 





0063 


0.0445 ± 


0010 





3412 


± 





0042 


63 ± 10 


19 ±3 


636 ± 1 


NTH45PH30 





0457 


± 





0001 


-0.0334 


± 





0066 


0.0395 ± 


0021 





3419 


± 





0041 


40 ± 10 


20 ±3 


820 ± 1 


NTH45PH40 





0459 


± 





0001 


-0.0397 


± 





0064 


0.0331 ± 


0031 





3422 


± 





0045 


20 ±8 


16 ±3 


967 ± 1 


NTH45PH45 





0459 


± 





0001 


-0.0425 


± 





0063 


0.0297 ±0 


0035 





3425 


± 





0046 


12 ±8 


12 ±3 


1026 ± 1 


NTH45PH50 





0460 


± 





0001 


-0.0455 


± 





0068 


0.0261 ± 


0038 





3426 


± 





0049 


3±8 


5±3 


1089 ± 1 


NTH45PH60 





0462 


± 





0001 


-0.0491 


± 





0062 


0.0177 ±0 


0049 





3432 


± 





0051 


-9 ±7 


-10±4 


1191 ± 2 


NTH45PH75 





0464 


± 





0001 


-0.0535 


± 





0069 


0.0056 ± 


0048 





3445 


± 





0049 


-18 ±7 


-41 ±4 


1318 ±3 


NTH45PH90 





0466 


± 





0001 


-0.0548 


± 





0074 


-0.0081 ± 


0.0061 





3455 


± 





0045 


-11±8 


-74 ±4 


1380 ± 3 


NTH45PH975 





0466 


± 





0001 


-0.0530 


± 





0068 


-0.0150 ± 


0.0065 





3453 


± 





0048 


-1±7 


-89 ±6 


1378 ± 3 


NTH45PH120 





0466 


± 





0001 


-0.0441 


± 





0052 


-0.0338 ± 


0.0075 





3452 


± 





0043 


47 ±9 


-111 ±6 


1217± 3 


NTH45PH1125 





0466 


± 





0001 


-0.0480 


± 





0059 


-0.0281 ± 


0.0075 





3454 


± 





0045 


29 ±9 


-108 ±6 


1297 ± 3 


NTH45PH1275 





0465 


± 





0001 


-0.0394 


± 





0042 


-0.0388 ± 


0.0073 





3448 


± 





0042 


65 ± 10 


-110±7 


1114± 3 


NTH45PH150 





0460 


± 





0001 


-0.0211 


± 





0003 


-0.0494 ± 


0.0058 





3431 


± 





0037 


110 ± 12 


-76 ±6 


633 ±2 


NTH60PH0 





0440 


± 





0001 


-0.0101 


± 





0070 


0.0640 ± 


0008 





3260 


± 





0042 


-109 ±4 


121 ± 1 


1435 ± 5 


NTH60PH30 





0444 


± 





0001 


-0.0449 


± 





0119 


0.0525 ± 


0003 





3271 


± 





0048 


-202 ± 2 


46 ± 1 


1568 ±4 


NTH60PH60 





0443 


± 





0001 


-0.0674 


± 





0138 


0.0269 ± 


0015 





3275 


± 





0044 


-194 ± 3 


-71 ±3 


1447 ± 3 


NTH60PH90 





0438 


± 





0000 


-0.0736 


± 





0141 


-0.0051 ± 


0.0015 





3263 


± 





0036 


-97 ±3 


-145 ±3 


1050 ± 2 


NTH60PH120 





0430 


± 





0000 


-0.0563 


± 





0063 


-0.0385 ± 


0.0044 





3233 


± 





0035 


25 ±4 


-96 ±3 


227 ± 1 


NTH60PH150 





0430 


± 





0000 


-0.0260 


± 





0010 


-0.0600 ± 


0.0041 





3231 


± 





0039 


26 ±4 


43 ±3 


-743 ± 3 


NTH120PH0 





0333 


± 





0000 


0.0040 ± 





0049 


0.0713 ±0 


0008 





2683 


± 





0012 


-105 ± 1 


48 ±2 


480 ±3 


NTH120PH30 





0330 


± 





0000 


-0.0339 


± 





0062 


0.0659 ± 


0018 





2667 


± 





0013 


-43 ±2 


-117 ± 3 


-101 ± 1 


NTH120PH60 





0332 


± 





0000 


-0.0617 


± 





0040 


0.0433 ± 


0015 





2678 


± 





0014 


108 ±3 


-145 ±4 


-516 ±3 


NTH120PH90 





0335 


± 





0001 


-0.0748 


± 





0021 


0.0085 ± 


0009 





2710 


± 





0005 


229 ±4 


-40 ±4 


-835 ± 5 


NTH120PH120 





0340 


± 





0001 


-0.0668 


± 





0007 


-0.0289 ± 


0.0005 





2736 


± 





0001 


223 ± 3 


147 ±3 


-1155 ±5 


NTH120PH150 





0339 


± 





0001 


-0.0401 


± 





0035 


-0.0589 ± 


0.0010 





2719 


± 





0006 


26 ±2 


199 ±2 


-1008 ±5 


NTH135PH0 





0315 


± 





0000 


0.0120 ± 





0016 


0.0597 ±0 


0008 





2576 


± 





0016 


-103 ± 1 


-105 ±2 


-670 ± 3 


NTH135PH30 





0318 


± 





0000 


-0.0204 


± 





0033 


0.0571 ± 


0001 





2582 


± 





0021 


21 ± 2 


-122 ±2 


-898 ± 5 


NTH135PH60 





0318 


± 





0000 


-0.0464 


± 





0031 


0.0389 ± 


0009 





2582 


± 





0025 


95 ±3 


-36 ± 1 


-915 ±6 


NTH135PH90 





0317 


± 





0000 


-0.0594 


± 





0017 


0.0111 ± 


0005 





2575 


± 





0031 


66 ±3 


75 ±1 


-772 ± 6 


NTH135PH120 





0314 


± 





0000 


-0.0570 


± 





0008 


-0.0197 ± 


0.0004 





2562 


± 





0025 


-61 ± 1 


120 ±2 


-370 ± 4 


NTH135PH150 





0313 


± 





0000 


-0.0402 


± 





0001 


-0.0455 ± 


0.0008 





2560 


± 





0019 


-154 ± 1 


22 ± 1 


210 ± 1 


NTH165PH0 





0292 


± 





0001 


0.0083 ± 





0017 


0.0220 ± 


0008 





2525 


± 





0023 


-98 ± 1 


142 ±3 


-288 ± 2 


NTH165PH30 





0292 


± 





0001 


-0.0044 


± 





0007 


0.0212 ±0 


0005 





2527 


± 





0022 


-82 ±1 


141 ±3 


-333 ± 2 


NTH165PH60 





0292 


± 





0001 


-0.0158 


± 





0006 


0.0146 ± 


0016 





2527 


± 





0022 


-75 ±1 


155 ±3 


-293 ± 2 


NTH165PH90 





0292 


± 





0001 


-0.0230 


± 





0016 


0.0041 ± 


0025 





2526 


± 





0021 


-84 ± 1 


170 ±3 


-176 ± 1 


NTH165PH120 





0292 


± 





0001 


-0.0240 


± 





0023 


-0.0074 ± 


0.0025 





2523 


± 





0023 


-100 ± 1 


171 ±3 


-8± 1 


NTH165PH150 





0292 


± 





0001 


-0.0187 


± 





0023 


-0.0169 ± 


0.0019 





2523 


± 





0024 


-107 ± 1 


157 ±3 


164 ± 1 
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TABLE XV: Radiated mass, angular momentum, and the remnant recoil (in original frame) as calculated from 1/14,. Errors 
quoted are from differences between to extrapolation to r = oo. See Tables [XVI| and |XVII| for more accurate measurement of 
the error. 



CONF 


5M rad 






SJx 






SJy 


SJ Z 






V x 


Vy 


V, 




STH45PH0 


0.0380 


± 


0.0001 


-0.0039 


± 


0.0022 


0.0033 ± 0.0029 


0.3026 


± 


0.0035 


-61 ± 1 


255 ±3 


783 ±3 




STH45PH30 


0.0371 


± 


0.0001 


-0.0046 


± 


0.0034 


0.0018 ± 0.0009 


0.2988 


± 


0.0036 


-14 ± 1 


261 ±2 


-471 ± 3 




STH45PH60 


0.0371 


± 


0.0001 


-0.0047 


± 


0.0043 


0.0013 ± 0.0008 


0.3000 


± 


0.0032 


-16 ± 2 


265 ±2 


-1335 ± 


6 


STH45PH90 


0.0376 


± 


0.0001 


-0.0035 


± 


0.0040 


0.0006 ±0.0010 


0.3036 


± 


0.0025 


-47 ±3 


265 ± 1 


-1835 ± 


8 


STH45PH120 


0.0384 


± 


0.0001 


-0.0003 


± 


0.0017 


-0.0016 ±0.0011 


0.3065 


± 


0.0025 


-70 ±3 


262 ± 1 


-2004 ± 


8 


STH45PH150 


0.0386 


± 


0.0001 


0.0028 ± 





0009 


-0.0032 ± 0.0030 


0.3068 


± 


0.0030 


-89 ±1 


254 ±2 


-1707 ± 


5 


KTH45PH0 


0.0389 


± 


0.0001 


-0.0026 


± 


0.0029 


0.1083 ±0.0020 


0.2797 


± 


0.0025 


147 ±2 


-486 ± 1 


-1777 ± 


4 


KTH45PH30 


0.0392 


± 


0.0001 


-0.0670 


± 


0.0145 


0.0932 ± 0.0006 


0.2802 


± 


0.0031 


477 ±2 


-314 ±2 


-2022 ± 


3 


KTH45PH60 


0.0394 


± 


0.0001 


-0.1143 


± 


0.0243 


0.0517 ±0.0020 


0.2809 


± 


0.0043 


578 ±4 


-7± 1 


-2042 ± 


5 


KTH45PH90 


0.0401 


± 


0.0001 


-0.1251 


± 


0.0207 


-0.0043 ± 0.0024 


0.2841 


± 


0.0053 


469 ±7 


327 ±2 


-1991 ± 


3 


KTH45PH105 


0.0408 


± 


0.0001 


-0.1166 


± 


0.0157 


-0.0336 ± 0.0030 


0.2873 


± 


0.0050 


295 ±5 


447 ±4 


-1805 ± 


2 


KTH45PH120 


0.0413 


± 


0.0001 


-0.1011 


± 


0.0115 


-0.0591 ± 0.0029 


0.2892 


± 


0.0048 


63 ±3 


430 ±5 


-1355 ± 


2 


KTH45PH135 


0.0412 


± 


0.0001 


-0.0809 


± 


0.0087 


-0.0795 ± 0.0028 


0.2884 


± 


0.0044 


-134 ± 1 


218 ±3 


-578 ± 1 




KTH45PH150 


0.0402 


± 


0.0001 


-0.0565 


± 


0.0059 


-0.0961 ± 0.0045 


0.2835 


± 


0.0047 


-186 ± 1 


-113 ±2 


411 ±3 




KTH45PH165 


0.0392 


± 


0.0001 


-0.0290 


± 


0.0024 


-0.1062 ± 0.0045 


0.2802 


± 


0.0036 


-64 ± 1 


-387 ± 1 


1280 ± 4 




KTH22.5PH0 


0.0385 


± 


0.0001 


-0.0025 


± 


0.0038 


0.0642 ± 0.0020 


0.2929 


± 


0.0038 


266 ±9 


-43 ± 1 


-1641 ± 


4 


KTH22.5PH30 


0.0387 


± 


0.0001 


-0.0320 


± 


0.0026 


0.0537 ±0.0015 


0.2943 


± 


0.0031 


337 ±9 


130 ± 1 


-1315 ± 


6 


KTH22.5PH60 


0.0386 


± 


0.0001 


-0.0531 


± 


0.0012 


0.0285 ± 0.0047 


0.2943 


± 


0.0028 


231 ±9 


296 ± 1 


-793 ± 5 




KTH22.5PH90 


0.0380 


± 


0.0001 


-0.0604 


± 


0.0001 


-0.0030 ± 0.0058 


0.2920 


± 


0.0034 


33 ±8 


315 ±1 


-23 ±3 




KTH22.5PH120 


0.0375 


± 


0.0001 


-0.0511 


± 


0.0017 


-0.0348 ± 0.0070 


0.2898 


± 


0.0037 


-71 ±10 


136 ± 1 


948 ± 1 




KTH22.5PH150 


0.0378 


± 


0.0001 


-0.0278 


± 


0.0039 


-0.0567 ± 0.0050 


0.2903 


± 


0.0044 


58 ±9 


-55 ±1 


1594 ± 2 




LPHO 


0.0538 


± 


0.0002 


-0.0085 


± 


0.0084 


0.0676 ± 0.0004 


0.3813 


± 


0.0071 


74 ± 10 


334 ±6 


2734 ± 3 




LPH30 


0.0544 


± 


0.0002 


-0.0460 


± 


0.0141 


0.0566 ±0.0011 


0.3834 


± 


0.0068 


-30 ±7 


246 ±5 


2181 ± 1 




LPH60 


0.0545 


± 


0.0002 


-0.0688 


± 


0.0132 


0.0312 ±0.0024 


0.3840 


± 


0.0064 


-51 ±9 


150 ±5 


1479 ± 1 




LPH90 


0.0538 


± 


0.0002 


-0.0765 


± 


0.0123 


-0.0050 ± 0.0053 


0.3819 


± 


0.0061 


36 ± 14 


95 ±5 


311 ±3 




LPH120 


0.0524 


± 


0.0002 


-0.0639 


± 


0.0083 


-0.0371 ± 0.0033 


0.3767 


± 


0.0068 


207 ± 14 


169 ±5 


-1473 ± 


3 


LPH150 


0.0525 


± 


0.0002 


-0.0327 


± 


0.0004 


-0.0614 ± 0.0021 


0.3771 


± 


0.0070 


231 ± 13 


341 ±6 


-2835 ± 


1 
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TABLE XVI: Comparison between remnant horizon properties and radiated quantities. Differences between the two is a much 
better measurement of the true error. 



CONF 


<5M rad 


5M rcm 




<5 J xrcm 


$J yrad 


S J yrcm 


3 J zrad 


$J zrem 




NTH15PH0 





0494 





0510 


-0.0015 


0.0059 


o 


0199 





0301 


o 


3690 


o 


3774 




NTH15PH30 





0493 





0509 


-0.0090 


-0.0178 


o 


0154 





0132 


o 


3687 


o 


3774 




NTH15PH60 


o 


0492 


o 


0508 


-0.0140 


-0.0266 


o 


0068 


o 


0065 


o 


3685 


o 


3770 




NTH15PH90 


o 


0492 


o 


0508 


-0.0153 


-0.0284 




0.0035 




0.0037 


o 


3690 


o 


3770 




NTH15PH120 


o 


0493 





0509 


-0.0126 


-0.0228 




0.0129 




0.0142 


o 


3689 


o 


3773 




NTH15PH150 


o 


0494 


o 


0510 


-0.0065 


0.0006 




0.0189 




0.0135 


o 


3692 


o 


3783 




NTH30PH0 


o 


0481 





0496 


-0.0040 


-0.0103 


o 


0377 





0419 


o 


3588 


o 


3673 




NTH30PH30 


o 


0478 





0493 


-0.0194 


-0.0079 


o 


0285 





0276 


o 


3576 


o 


3676 




NTH30PH60 


o 


0477 





0493 


-0.0288 


-0.0419 


o 


0120 





0186 


o 


3575 


o 


3661 




NTH30PH90 


o 


0480 





0495 


-0.0324 


-0.0426 




0.0075 




0.0047 


o 


3587 


o 


3674 




NTH30PH120 


o 


0483 


o 


0499 


-0.0271 


-0.0207 




0.0252 




0.0262 


o 


3602 


o 


3697 




NTH30PH150 


o 


0483 





0499 


-0.0140 


-0.0234 




0.0362 




0.0322 


o 


3601 


o 


3687 




NTH45PH0 


o 


0456 





0470 


-0.0075 


-0.0052 


o 


0507 





0588 


o 


3410 


o 


3496 




NTH45PH10 


o 


0456 





0470 


-0.0167 


-0.0183 


o 


0481 





0399 


o 


3409 


o 


3520 




NTH45PH20 


o 


0456 


o 


0471 


-0.0258 


-0.0215 


o 


0445 





0380 


o 


3412 


o 


3523 




NTH45PH30 


o 


0457 





0472 


-0.0334 


-0.0329 


o 


0395 





0528 


o 


3419 


o 


3504 




NTH45PH40 


o 


0459 





0473 


-0.0397 


-0.0358 


o 


0331 





0463 


o 


3422 


o 


3515 




NTH45PH45 


o 


0459 


o 


0474 


-0.0425 


-0.0413 


o 


0297 


o 


0420 


o 


3425 


o 


3518 




NTH45PH50 


o 


0460 





0475 


-0.0455 


-0.0488 


o 


0261 





0398 


o 


3426 


o 


3517 




NTH45PH60 


o 


0462 





0476 


-0.0491 


-0.0436 


o 


0177 





0089 


o 


3432 


o 


3552 




NTH45PH75 


o 


0464 





0479 


-0.0535 


-0.0563 


o 


0056 





0187 


o 


3445 





3538 




NTH45PH90 


o 


0466 





0481 


-0.0548 


-0.0523 




0.0081 




0.0145 


o 


3455 


o 


3556 




NTH45PH975 


o 


0466 


o 


0481 


-0.0530 


-0.0497 




0.0150 




0.0188 


o 


3453 


o 


3557 




NTH45PH120 


o 


0466 





0481 


-0.0441 


-0.0393 




0.0338 




0.0413 


o 


3452 





3547 




NTH45PH1125 


o 


0466 


o 


0481 


-0.0480 


-0.0455 




0.0281 




0.0317 


o 


3454 


o 


3551 




NTH45PH1275 


o 


0465 





0480 


-0.0394 


-0.0337 




0.0388 




0.0458 


o 


3448 


o 


3542 




NTH45PH150 


o 


0460 





0475 


-0.0211 


-0.0215 




0.0494 




0.0548 


o 


3431 


o 


3514 




NTH60PH0 


o 


0440 





0455 


-0.0101 


-0.0120 


o 


0640 





0716 


o 


3260 


o 


3347 




NTH60PH30 


o 


0444 


o 


0459 


-0.0449 


-0.0279 


o 


0525 


o 


0521 


o 


3271 


o 


3384 




NTH60PH60 


o 


0443 





0457 


-0.0674 


-0.0522 


o 


0269 





0308 


o 


3275 


o 


3378 




NTH60PH90 


o 


0438 





0452 


-0.0736 


-0.0521 




0.0051 





0000 


o 


3263 


o 


3372 




NTH60PH120 


o 


0430 





0443 


-0.0563 


-0.0613 




0.0385 




0.0366 


o 


3233 


o 


3305 




NTH60PH150 


o 


0430 


o 


0444 


-0.0260 


-0.0362 




0.0600 




0.0527 


o 


3231 


o 


3319 




NTH120PH0 


o 


0333 


o 


0343 


0.0040 


0.0197 


o 


0713 





0722 


o 


2683 


o 


2732 




NTH120PH30 


o 


0330 


o 


0340 


-0.0339 


-0.0382 


o 


0659 


o 


0635 


o 


2667 


o 


2718 




NTH120PH60 


o 


0332 





0341 


-0.0617 


-0.0512 


o 


0433 





0453 


o 


2678 


o 


2745 




NTH120PH90 


o 


0335 





0344 


-0.0748 


-0.0675 


o 


0085 





0101 


o 


2710 


o 


2762 




NTH120PH120 


o 


0340 





0350 


-0.0668 


-0.0788 




0.0289 




0.0384 


o 


2736 


o 


2748 




NTH120PH150 


o 


0339 





0349 


-0.0401 


-0.0412 




0.0589 




0.0565 


o 


2719 


o 


2763 




NTH135PH0 


o 


0315 


o 


0325 


0.0120 


0.0112 


o 


0597 


o 


0560 


o 


2576 


o 


2634 




NTH135PH30 


o 


0318 





0327 


-0.0204 


-0.0187 


o 


0571 


o 


0564 


o 


2582 


o 


2646 




NTH135PH60 





0318 





0327 


-0.0464 


-0.0410 





0389 





0341 





2582 





2657 




NTH135PH90 





0317 





0326 


-0.0594 


-0.0615 





0111 





0206 





2575 





2639 




in inioornizu 





0314 





0323 


— u.uo / u 


— U.UDOO 




0.0197 




0.0092 





2562 





2622 




NTH135PH150 





0313 





0322 


-0.0402 


-0.0370 




0.0455 




0.0386 





2560 





2627 




NTH165PH0 





0292 





0301 


0.0083 


0.0143 





0220 





0291 





2525 





2576 




NTH165PH30 





0292 





0301 


-0.0044 


-0.0107 





0212 





0183 





2527 





2579 




NTH165PH60 





0292 





0301 


-0.0158 


-0.0211 





0146 





0116 





2527 





2577 




NTH165PH90 





0292 





0300 


-0.0230 


-0.0130 





0041 





0093 





2526 





2583 




NTH165PH120 





0292 





0300 


-0.0240 


-0.0176 




0.0074 





0046 





2523 





2583 




NTH165PH150 





0292 





0300 


-0.0187 


-0.0206 




0.0169 




0.0213 





2523 





2573 
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TABLE XVII: Comparison between remnant horizon properties and radiated quantities. Differences between the two is a much 
better measurement of the true error. 



CONF 


5M rad SM Icm 


foJxY&d $J a; rem 


^yrad $J yrcm 


$J zrad $J zrem 




STH45PH0 

STH45PH30 

STH45PH60 

STH45PH90 

STH45PH120 

STH45PH150 


0.0380 0.0396 
0.0371 0.0386 
0.0371 0.0386 
0.0376 0.0391 
0.0384 0.0400 
0.0386 0.0403 


-0.0039 0.0029 
-0.0046 0.0049 
-0.0047 0.0060 
-0.0035 -0.0044 
-0.0003 0.0021 
0.0028 0.0007 


0.0033 0.0113 
0.0018 0.0106 
0.0013 0.0059 
0.0006 -0.0119 
-0.0016 0.0083 
-0.0032 -0.0130 


0.3026 0.3100 
0.2988 0.3066 
0.3000 0.3071 
0.3036 0.3097 
0.3065 0.3127 
0.3068 0.3134 




KTH45PH0 

KTH45PH30 

KTH45PH60 

KTH45PH90 

KTH45PH105 

KTH45PH120 

KTH45PH135 

KTH45PH150 

KTH45PH165 


0.0389 0.0407 
0.0392 0.0410 
0.0394 0.0410 
0.0401 0.0417 
0.0408 0.0424 
0.0413 0.0430 
0.0412 0.0429 
0.0402 0.0420 
0.0392 0.0410 


-0.0026 -0.0097 
-0.0670 -0.0624 
-0.1143 -0.0927 
-0.1251 -0.1215 
-0.1166 -0.1208 
-0.1011 -0.0903 
-0.0809 -0.0688 
-0.0565 -0.0628 
-0.0290 -0.0360 


0.1083 0.1226 
0.0932 0.1065 
0.0517 0.0556 
-0.0043 0.0020 
-0.0336 -0.0407 
-0.0591 -0.0560 
-0.0795 -0.0758 
-0.0961 -0.0995 
-0.1062 -0.1073 


0.2797 0.2830 
0.2802 0.2851 
0.2809 0.2912 
0.2841 0.2903 
0.2873 0.2915 
0.2892 0.2990 
0.2884 0.2990 
0.2835 0.2901 
0.2802 0.2869 




KTH22.5PH0 

KTH22.5PH30 

KTH22.5PH60 

KTH22.5PH90 

KTH22.5PH120 

KTH22.5PH150 


0.0385 0.0401 
0.0387 0.0403 
0.0386 0.0402 
0.0380 0.0396 
0.0375 0.0391 
0.0378 0.0394 


-0.0025 0.0103 
-0.0320 -0.0232 
-0.0531 -0.0497 
-0.0604 -0.0621 
-0.0511 -0.0589 
-0.0278 -0.0388 


0.0642 0.0695 
0.0537 0.0594 
0.0285 0.0387 
-0.0030 0.0085 
-0.0348 -0.0412 
-0.0567 -0.0611 


0.2929 0.3018 
0.2943 0.3036 
0.2943 0.3030 
0.2920 0.3010 
0.2898 0.2976 
0.2903 0.2988 




LPHO 

LPH30 

LPH60 

LPH90 

LPH120 

LPH150 


0.0538 0.0562 
0.0544 0.0568 
0.0545 0.0567 
0.0538 0.0560 
0.0524 0.0546 
0.0525 0.0548 


-0.0085 0.0001 
-0.0460 -0.0336 
-0.0688 -0.0657 
-0.0765 -0.0650 
-0.0639 -0.0494 
-0.0327 -0.0280 


0.0676 0.0788 
0.0566 0.0663 
0.0312 0.0233 
-0.0050 0.0082 
-0.0371 -0.0439 
-0.0614 -0.0688 


0.3813 0.3941 
0.3834 0.3969 
0.3840 0.3977 
0.3819 0.3956 
0.3767 0.3906 
0.3771 0.3902 
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TABLE XVIII: BH spins during final plunge, recoil velocity, and the angle between Ax for PHYYY and Ax of the corresponding 
PHO configuration; all calculated in a rotated frame where the infall occurs in the xy plane. 



CONF 


S x l 


S 


yl 


S z l 


S 


x2 


Sy2 


s 


z2 


14 


Vy 


v z 


V 


NTH15PH0 


0.0319 
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0.2000 





0000 





0000 





0000 


77 


-110 


-457 





NTH15PH30 
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0.0238 


0.1998 
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0000 





0000 


79 
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-252 


24.5 


NTH15PH60 


0.0470 




0.0014 


0.2012 





0000 





0000 





0000 


82 


-117 


15 


50.5 


NTH15PH90 


0.0372 





0210 


0.2022 





0000 





0000 





0000 


85 


-126 


279 


81.7 


NTH15PH120 


-0.0249 




0.0216 


0.2030 





0000 





0000 





0000 


-116 


-98 


476 


122.8 
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79 


-117 
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84 
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0000 


66 
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0000 
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47.4 
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0.0645 
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0.1646 





0000 





0000 





0000 
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-109 
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56.3 
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0.0368 
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0.1615 





0000 





0000 





0000 


100 


-93 


1311 


70.9 
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0.0061 





1328 
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0000 





0000 





0000 


82 


-72 


1377 


85.1 
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-0.0093 





1353 


0.1552 





0000 





0000 





0000 


72 


-62 


1377 


91.6 
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-0.0640 





1261 


0.1498 





0000 





0000 





0000 


47 


-29 


1222 


114.6 


NTH45PH1125 


-0.0430 





1331 


0.1512 





0000 





0000 





0000 


52 


-41 


1300 


105.6 


NTH45PH1275 


-0.0687 





1234 


0.1497 





0000 





0000 





0000 


38 


-37 


1120 


116.8 


NTH45PH150 


-0.1109 





0837 


0.1516 





0000 





0000 





0000 


46 


-44 


644 


140.7 


NTH60PH0 


0.0666 





1442 


0.1310 





0000 





0000 





0000 


103 


-72 


1439 





NTH60PH30 


-0.0104 





1677 


0.1188 





0000 





0000 





0000 


65 


-11 


1581 


28.3 


NTH60PH60 


-0.0465 





1658 


0.1127 





0000 





0000 





0000 


35 


6 


1461 


40.4 


NTH60PH90 


-0.1169 





1294 


0.1088 





0000 





0000 





0000 


17 


33 


1064 


66.9 


NTH60PH120 


-0.1616 





0407 


0.1195 





0000 





0000 





0000 


54 


-15 


242 


100.7 


NTH60PH150 


-0.1399 




0.0672 


0.1345 





0000 





0000 





0000 


113 


-85 


-731 


140.4 


NTH120PH0 


-0.1641 





0877 


-0.0856 





0000 





0000 





0000 


7 


152 


470 





NTH120PH30 


-0.1909 




0.0090 


-0.0731 





0000 





0000 





0000 


27 


104 


-119 


30.8 


NTH120PH60 


-0.1803 




0.0696 


-0.0682 





0000 





0000 





0000 


47 


77 


-539 


49.2 


NTH120PH90 


-0.1519 




0.1190 


-0.0694 





0000 





0000 





0000 


54 


71 


-862 


66.2 


NTH120PH120 


-0.0569 




0.1801 


-0.0833 





0000 





0000 





0000 


38 


113 


-1179 


100.6 


NTH120PH150 


0.0757 




0.1666 


-0.0944 





0000 





0000 





0000 


4 


168 


-1014 


142.6 


NTH135PH0 


-0.1314 




0.1005 


-0.1230 





0000 





0000 





0000 


45 


122 


-673 





NTH135PH30 


-0.0650 




0.1461 


-0.1307 





0000 





0000 





0000 


37 


143 


-894 


28.6 


NTH135PH60 


0.0019 




0.1537 


-0.1371 
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o 


0000 





0000 


25 


169 


-905 


53.3 


NTH135PH90 


0.0664 




0.1355 


-0.1383 





0000 





0000 





0000 


15 


184 


-757 


78.7 


NTH135PH120 


0.1391 




0.0699 


-0.1354 





0000 





0000 





0000 


12 


174 


-352 


115.9 


NTH135PH150 


0.1620 





0274 


-0.1237 





0000 





0000 





0000 


30 


138 


221 


152.2 


NTH165PH0 


-0.0416 




0.0470 


-0.1949 





0000 





0000 





0000 


33 


194 


-272 





NTH165PH30 


-0.0132 




0.0579 


-0.1992 





0000 





0000 





0000 


32 


197 


-313 


28.7 


NTH165PH60 


0.0197 




0.0533 


-0.1999 





0000 





0000 





0000 


29 


200 


-273 


61.8 


NTH165PH90 


0.0419 




0.0385 


-0.1993 





0000 





0000 





0000 


23 


200 


-162 


89. 


NTH165PH120 


0.0625 




0.0059 


-0.1992 





0000 





0000 





0000 


28 


197 


-3 


126.2 


NTH165PH150 


-0.1161 





0870 


-0.1459 





0000 





0000 





0000 


-87 


-126 


199 


153.2 
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TABLE XIX: BH spins during final plunge, recoil velocity, and the angle between Ai for PHYYY and Ax of the corresponding 
PHO configuration; all calculated in a rotated frame where the infall occurs in the xy plane. 
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TABLE XX: Comparing the measured value of the z component of the recoil to values of Ax , S± , A z , and S z in a rotated 
frame where the infall occurs in the xy plane. These data are equivalent to the data given in Table |XiX"j 
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